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Introduction 


I offer this work as the mathematical principles of philosophy, for the whole 
burden of philosophy seems to consist in this — from the phenomena of 
motions to investigate the forces of nature, and then from these forces to 
demonstrate the other phenomena; and to this end the general propositions in 
the first and second Books are directed. In the third Book 1... derive from 
celestial phenomena the forces of gravity with which bodies tend to the Sun 
and the several planets. Then from these forces, by other propositions which 
are also mathematical, I deduce the motions of the planets, the comets, the 
Moon, and the sea. 


Sir Isaac Newton, from the Preface to the First Edition of Principia. 


The third book of Sir Isaac Newton’s monumental work Principia, entitled The System 
of the World, was published in 1687, some twenty years after he had discovered many of 
its results. In it he establishes his universal law of gravitation and uses his three laws of 
motion to predict the character of the orbits of planets and their moons. Tradition has 
it that Halley, the astronomer, Hooke, the physicist, and Wren, the architect, had 
discussed at various times in 1684 the problem of determining the orbit of a planet 
around the Sun if the force between them varied as the inverse square of the distance. 
Unable to succeed with the problem, they took it to Newton, only to find that he had 
already solved it, together with other problems of planetary motion. Encouraged 
constantly by Halley, Newton finally published Principia. The theory which it 
expounded was extraordinarily successful. Only in the twentieth century, with the 
theories of relativity and quantum mechanics, has it become evident that there are 
natural processes not explained accurately by the Newtonian world-view. 


The main aim of this unit is to apply Newtonian mechanics in order to derive Kepler’s 
three laws of planetary orbits. Johann Kepler (1571-1630) arrived at these laws by trial 
and error, from careful analysis of the measurements of planetary orbits taken by 
himself and by Tycho Brahe (1546-1601), the Danish astronomer. The first two laws 
were published in 1609 and the third in 1619, so they precede Newton’s mathematical 
derivations by a good lifetime. 


In deriving the laws of planetary motion we shall use modern notation and methods, 
although it must be said that the calculus in its original form was the innovation of 
Newton and of his contemporary, Leibnitz. The modern use of geometric vectors was 
apparently an invention of J. Willard Gibbs (1839-1903), and is a major notational 
advance over the geometric methods which held sway in Newton’s era. 


In summary, then, there is little new mathematics introduced in this unit. Instead we 
utilize the methods of calculus and vector analysis developed in the course to describe 
one of the great scientific advances of the last three centuries. 


Study guide 


The material in Section 1 is largely revision, involving methods which will be of 
particular use later in the unit. A review of energy conservation for one-dimensional 
motion (from Unit 7) leads on to the corresponding topic for a particle moving in three 
dimensions. Then the planar motion of a particle is described in terms of polar 
coordinates and the associated unit vectors. This is a generalization of the treatment in 
Unit 28, where the particle motion was confined to a circle. 


The time spent on Section 1 will depend upon your confidence with the previous 
- mechanics in the course, but you should not spend more than 30% of your total study 
time here. The major part of your effort should be reserved for Sections 2 and 3. 


Kepler’s laws for planetary orbits are stated and explained in Section 2. Following this, 
Newton’s universal law of gravitation is introduced. In Section 3 the equation of 
gravitational orbits is derived, and it is shown how Kepler’s laws follow as a 
consequence of Newtonian mechanics. 
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The television programme in Section 4 covers a mathematical development similar to 
that in Section 3, but with greater reference to non-gravitational examples. More 
benefit will probably be obtained from the programme if you have managed to complete 
most of Section 3 first, but it should still be of help if viewed at an earlier stage. 


1 Review of energy and planar motion 


The ultimate goal of this unit is to model the motions of planets, satellites and comets. 
To facilitate this study it is useful to review some aspects of Newtonian mechanics, and 
that is the purpose of Section 1. 


Subsection 1.1 will remind you about the material in Unit 7 concerning particle 
motions in one dimension for which the total mechanical energy is constant. In 
Subsection 1.2 the corresponding motions in three dimensions are considered, and it is 
shown that the total mechanical energy is conserved provided that the force acting on 
the particle is conservative in the sense of Unit 26. This conclusion applies also in one 
or two dimensions. 


Subsection 1.3 shows how motion in the plane may be described in terms of polar 
coordinates r,@ and the associated unit vectors e,,eg. You saw circular motion 
expressed in this way in Unit 28. Following the derivation of kinematical formulas, 
there are some applications of Newton’s second law to find equations of motion in this 
coordinate system. 


1.1 Mechanical energy in one dimension 


Consider a particle of mass m, moving in one dimension under the action of a force. We 
choose the z-axis to coincide with the line along which the particle moves. Then, as 
you know from Unit 7, Newton’s second law is mi = F’, where F is the x-component of 
the force. In general, the force may depend upon the location x of the particle, its 
velocity z, the time t, or any combination of these variables. We consider here the case 
where F’ depends only upon the particle’s position x. (A special case of this is when F' 
is constant.) A position-dependent force in one dimension is called conservative 
because, as you may recall from Unit 7 and as we shall shortly verify anew, the total 
mechanical energy is constant (conserved) for such a force. 


In terms of the force F' = F(x), the potential energy U = U(x) of the particle is defined 
by the equation 


= — (1) 


where Zo is any conveniently chosen point. 


Exercise 1 


(i) Find the potential energy U(x) corresponding to the force F(x) = k, where k is a 
constant. 


(ii) Find the potential energy U(x) corresponding to the force F(a) = —kx, where k is a 
constant. 


(iii) 


Find the force F(x) corresponding to the potential energy U(x) = k/zx (x > 0), where k is 
a constant. 


(Solution on page 47] 


When it is possible to define a potential energy function then the sum of the kinetic 
energy smi? and the potential energy U(x) is constant throughout the motion of the 
particle. The sum : 


E = imz* + U(z) | (2) 


Unit 7 Section 3 


Unit 26 Subsection 4.3 


Unit 28 Subsection 4.2 


You may see a particle 
described elsewhere as either 
a ‘point mass’ or a ‘mass 
point’. 


The description 
‘conservative’ was not used 
when energy was introduced 
in Unit 7 Section 3. You will 
see in the next subsection 
that energy-conserving forces 
are also conservative in the 
sense of Unit 26. 


The point xo is the datum of 
the potential energy function 
U, for which U(2o) = 0. This 
point (or equivalently, the 
constant of integration 
arising from the indefinite 
integral for U) is usually 
chosen so as to achieve the 
simplest possible form for 
lite). 


6 MST204 30.1 


is called the total mechanical energy of the particle, and its constancy is expressed by 
the equation F = 0. To establish that E is constant, we differentiate Equation (2) with 
respect to time ¢t. Using the chain rule on the right-hand side, this leads to 


= eae Hy : *) 
E=mzaz + —r2=2z(|mz+ —  ]. 
dx dx 


But the term in brackets vanishes since, by Newton’s second law and Equation (1), 


Hence we have E = 0. The fact that E is constant for a conservative force often 
provides useful information about the particle motion. For example, consider the 
harmonic oscillator (that is, a particle attached to one end of a perfect spring whose 
other end is fixed). Measuring x from the equilibrium position of the particle, the force 
F is given by , 


where k is the stiffness of the spring. The potential energy function is therefore 


U(x) = Lkx® + C, 


moe 
where C’ is a constant. Any value of C may be chosen here since the value of 
F(x) = —dU/dz is unaffected by the choice. One often takes C = 0. Figures 1(a) 
and 1(b) below show the potential energy U(x) for a negative constant C and for the 
more conventional choice C' = 0 respectively. 


U Bie 


&: 
(a) (b) 
Figure 1 


Either form for U(x) gives the correct equation of motion. For any chosen value of C 
the total mechanical energy EF = mz? + U(x) is a constant, although its value 
depends on the choice of C. 


To see how useful this knowledge is, suppose that initially, at t = 0, the particle has 
position Zp and is moving with velocity vp. This information establishes the total 
mechanical energy of the particle for all later times t, since E(t) = E(0). Whatever 
value of C' is chosen, we have 


E = mvp + 4kag+C = sme" (t) + 4kx?(t) +C. 


Solving for the speed squared gives 
2 2 
p(t) = — (E — 1ka?(t)-C) = —(B- 
i (t) = — (E - 3he*(t) -C) = = (B-U(2) 


Since the left-hand side can never be negative, we know that the term in brackets on 
the right-hand side is also non-negative. Hence, for all times t we know that 


Ulin) = 8. 


This information is conveyed graphically in Figure 2, where we have chosen the value of 
C’ used in Figure 1(a). The constant value of E is represented by the horizontal dashed. 
line. At each value of x the vertical distance between this horizontal line and the curve 
for U(x) gives the kinetic energy $m? of the particle at that value of x. This 
difference, E — U(x), does not depend on C. When the particle is located at x = a or 
wt x = G3, its kinetic energy is zero. The particle approaches these points, comes Figure 2 
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momentarily to a halt, and then reverses its direction of motion. For this reason, points 
like a and G are called turning points. 


Note that the motion of the particle is periodic: at any point 7, between the turning 
points the speed is \/2(£ — U(x,))/m, and this value is always reproduced when the 
particle returns to that point. ‘The velocity at the point takes the values 

+,/2(E — U(z,))/m alternately, depending on whether the particle is passing through 
x, from left to right or from right to left. The particle passes through zx, twice in each 
complete cycle. 


In fact, you know from the explicit solution for simple harmonic motion in Unit 7 that Unit 7 Subsection 2.3 
in this case the particle oscillates sinusoidally between the turning points. 


Exercise 2 


The position function for a particle of mass m, undergoing simple harmonic motion due to a 
perfect spring of stiffness k, can be written as 


vo. 
x(t) = ro coswt + — sinwt, 
Ww 


where w = ,/k/m is the angular frequency and Zo, vo are respectively the position and velocity 
of the particle at time t = 0. Verify explicitly that the total mechanical energy 


B= ima” + tka +C 
is constant, and find the turning points for this particular motion. 


‘Solution on page 47 | 


Explicit formula solutions to Newton’s second law for the case of conservative forces 
cannot always be derived. If the potential energy function U(x) is known, however, 
then the approach used prior to Exercise 2 may be applied. This gives qualitative 
information about the particle motion even when the position function can be 
established only by numerical computation. In all cases it is useful to know that F is 
constant. Consider, for example, the motion of a particle under the influence of the 
potential energy function shown in Figure 3. 


Figure 3 
Given only the information provided by this graph, we can make the following 
statements about the particle’s motion if it has the energy F indicated in the figure. 


(i) The motion is periodic, and is confined between the turning points a and e. At 
these points the particle slows to a stop and reverses its direction of motion. 


(ii) At any point x for which a < x < «¢, the particle’s velocity is +,/2(k — U(x))/m. 


(iii) The particle is moving fastest at z = 3, since E — U (x) has an overall maximum 
there. 


(iv) As the particle passes through x = y in either direction, it slows down and then 
speeds up again. This follows from the fact that the speed ,/2(E — U(x))/m has a 
local minimum there. 


(v) As the particle passes through x = 6 it experiences a local maximum of speed. 
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Example 1 


What can be said about the motion of a particle under the influence of the potential 
energy function shown in Figure 4, if the energy £ of the particle is as indicated on the The potential energy 
figure? In particular, how many turning points are there at this energy level, what are _ function of Figure 4 is the 


these turning points and where does the particle move fastest? same as that in Figure 3, but 
the total mechanical energies 


of the particles are different. 


U 


Figure 4 


Solution 


The particle is confined (see Figure 5) to periodic motion either between points a and 
y or between points 6 and ¢, depending upon its initial position z(0). 


Figure 5 


For instance, if it were released from rest at point a then it would pursue a periodic 
motion between the turning points a and 4, and its maximum speed would occur at 

x = 3. Alternatively, it might be confined by its initial conditions to move between 6 
and @, with maximum speed at x = c. Motion at this energy F is not possible between 
points 7 and 6, as the ‘kinetic energy’ E — U (x) would be negative in this region. For 
the same reason, motion at the given energy is impossible for x <a orforz>d¢. O 


It is also possible with such energy diagrams to consider the motion of a particle which 
is not bound to periodic motion. Such unbound motion occurs for the case depicted in 
Figure 6, where U(x) = 0 for x > 6. 


For example, suppose that the particle were released from rest at x = a, which fixes its 
total mechanical energy for all subsequent times at E = 5mz*(0) +U(a) = U(a). The 
particle would then accelerate to the right until it picked up maximum speed at x = £. 
From z = ( to x = y it would decelerate to a local speed minimum at x = y. It would 
then accelerate again to the right until it reached x = 6 , after which it would continue 
to the right with constant speed, in unbound motion. The speed is constant for x > 6 
because the potential energy in that region is constant (zero, in fact). 
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Figure 6 


Exercise 3 


For the particle of Figure 6, released from rest at x = a, give in terms of m and U expressions 
for the particle’s speed at x = 3, x = y and for x > 6. Arrange the speeds at the points 3, 
and 6 in order of magnitude. 


Exercise 4 


A particle is projected to the left in Figure 6 from some initial point zo to the right of 6, and 
with the total mechanical energy EF shown in the figure. Give its initial velocity 7(0) in terms 
of FE, and describe its subsequent motion. 


Solutions on page 47 | 


1.2 Mechanical energy in three dimensions 


Just as for one-dimensional motion, the force acting on a particle moving in two or 
three dimensions can depend upon position r, velocity r and time t, or on any 
combination of these variables. It often happens, however, that the force F depends 
only upon position, which may be denoted by writing F = F(r) or F = F(z, y, z). If in 
addition F can be expressed as the gradient of a scalar field, then it is called a 
conservative force. By convention the scalar field function here is written as —U, 
where U = U(r) = U(z, y, z) is called the potential energy of the particle. Thus in 
three (or fewer) dimensions, a conservative force F and potential energy U are 
connected by the equation 


F = — gradU. (3) 


Note that in one dimension this reduces to 


which is Equation (1) of the previous subsection. 


Conservative forces were discussed briefly in Unit 26 Subsection 4.3. It is enough for 
the purposes of this unit to recognize that where the force acting on the particle can be 
written as in Equation (3), the total mechanical energy E is constant, where 


E = imr? + U(z,y, 2). (4) 


Here we have written r? as a shorthand for r-r, which is also equal to |r|?. In 


Equation (4) the kinetic energy smr?, the potential energy U and the total energy E 


are scalars. In order to show that E is constant, it suffices to establish that E = 0. 
Differentiation of Equation (4) with respect to time t and use of the chain rule produces 


sd See 
E= a (mr -¢+ U(a,y, z)) 
OU dx OU dy 


ee os Sere Se 
= Oe ay He 


=mr-r+r-gradU 
=r-(mr+gradU). 


aU dz 
Oz dt 


The gradient of a scalar field 
og, denoted by grad ¢ or V4, 
was introduced in Unit 26 
Section 2. It is defined as the 
vector field 
Od. Oo 
Oa Oy 
Conservative vector fields 
were defined in terms of line 
integrals in Unit 26 
Section 3. In Subsection 4.3 
of that unit it was shown 
that the definition is 
equivalent to the condition 
given here, involving the 
gradient of a scalar field. 


As pointed out in Unit 26 
Subsection 4.3, the potential 
energy U can be expressed as 
the line integral 


~ f Pear 


For the one-dimensional case 
this reduces to 


u=- | Fae, 


as in Subsection 1.1. 


The derivative of a scalar 
product was derived in 
Unit 14 Subsection 4.3. 


This form of the chain rule 
was introduced in Unit 25 
Subsection 2.4. 
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But by Equation (3), we have gradU = —F. Also F = m# by Newton’s second law, so 
the term in brackets vanishes and hence E = 0. We have shown that, when the force has 
the form given by Equation (3), the total mechanical energy is constant. This outcome 
represents the natural extension to three dimensions of the results in Subsection 1.1. 


Exercise 5 


A particle of mass m moves under the influence of a uniform gravitational field, so that 
F = —mgj, where the unit vector j is directed vertically upwards. Show that this field is 
conservative by finding the corresponding potential energy function U. 


Exercise 6 
Given that F is a conservative force with potential energy U = k /r, where k is a constant and 
r = (x? + y? + 27)1/?, show that 
Ze ; kr 
i= 73 (ti + yj + 2k) = 


[Solutions on page 47 | 


There is a useful generalization to the result above concerning the constancy of E. 
Suppose that the particle is being acted upon by a conservative force, as defined by 
Equation (3), plus any force G which is always directed perpendicular to the particle’s 
velocity. Then its equation of motion is 


mr = —gradU +G, 


where G-r = 0. In the exercise below you are asked to show that, under these more 
general circumstances, F is still constant. 


Exercise 7 


Suppose that a particle is acted upon by a net force — grad U + G, where G-r = 0. Show that 
if F is given by Equation (4) then EF = 0. 


[Solution on page 47 | 


The result of Exercise 7 is often applied when modelling the constrained motion of 
objects. Consider, for instance, a particle of mass m sliding without friction down an 
inclined plane (see Figure 7). Since we are assuming that there is no friction present, 
the only forces acting are gravity and the reaction force R. The latter is always 
perpendicular to the plane and hence also to the velocity. The gravitational force —mgj 
is conservative (as shown in Exercise 5) with potential energy function 


U =mgy + C, 


where C’ is any constant. If u is the speed then, from the result of Exercise 7, the total 


mechanical energy E is constant, where Figure 7 


E=imu?+mgy+C. 


You saw another example of 
this type, featuring a force 
Exercise 8 directed at right angles to 
the velocity, when studying 
the motion of a pendulum in 
Unit 15 Subsection 5.4. 


In the situation just described, the particle is released from rest at a vertical height h above 
the z-axis in Figure 7. Find its speed wo when it reaches the origin. 


[Solution on page 48 | 
1.3. Planar motion using polar coordinates 
In this subsection we consider a means of describing the motion of a particle in a plane. 


For some applications it is convenient to describe the position vector r of the particle in 
terms of Cartesian coordinates (x,y) and Cartesian unit vectors i, j, that is, 


r(t) = x(t)it y(d)j. (5) 
Then the velocity and acceleration vectors are respectively 

v(t) = #(t) = 2(t)i + g(t); (6) 
and a(t) = V(t) = £(t) = £()i+ g(t)j. (7) 


The ease of derivation of Equations (6) and (7) from Equation (5) is due to the fact 
that i and j are constant vectors throughout the motion of the particle. In summary, 
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then, the motion of a particle in a plane can be described by giving its Cartesian 
coordinates (x,y) as functions of time, together with the corresponding fixed 
orthogonal unit vectors i and Jj. 


In Unit 28 you saw that it is sometimes more convenient to use polar coordinates [r, 6] 
and the corresponding orthogonal unit polar vectors e, and eg (see Figure 3). For 
instance, corresponding to Equation (5), the polar form of the position vector is 


r= Te,, (8) 


where r = |r| is the distance of the particle from the origin. This equation is not as 
simple as it looks because, for general motion in the plane, both r and e, depend on 
time, as do the polar angle 6 and its associated unit vector eg. 


To see this time dependence, consider the fact that any vector in the plane can be 
expressed in terms of either the pair i,j or the pair e,,eg. In particular, each of i and j 
themselves can be expressed as a linear combination of e,,eg, and vice versa. From 
Unit 28, or by noting that the unit vectors e,,eg are obtained by rotating i, j 
respectively through an angle 6 (see Figure 9), these relations are 


e, = cos#i+sin6j, €g = —sin#i+cos6j (9) 


(10) 


Equations (9) show that e, and eg are functions of the polar angle 0. Thus, if a particle 
is moving along some curve so that @ is time-dependent, then e, and eg depend on time 
too. For instance, Figure 10 shows the curve followed by a particle moving anticlockwise 
around and away from the origin, which passes through the point [ry, 6,] at time t; and 
the point [r2, 42] at some later time t2. Note that, while e, and eg depend upon 6, they 
are independent of the radial distance r (which does not appear in Equations (9)). 


and i= cose, — sin# eg, j=sinde, + cos0eg. 


Figure 10 


Exercise 9 


Show explicitly from Equations (9) that, at each point other than the origin, e,,e9 form an 
orthogonal pair of unit vectors, that is, they satisfy the conditions 


e,-e9 =0 and e,-e, =e¢:-e, = l. 
(You will need to use the fact that i,j form an orthogonal pair of unit vectors.) 


[Solution on page 48 | 


In order to study the motion of a particle in plane polar coordinates, it is necessary to 
have expressions for its velocity and acceleration vectors in terms of r,6,e, and eg. In 
the course of deriving these equations, we need to express the time derivatives é,,é@9 in 
terms of e,,eg. The required relations for achieving this are 


e,. = Jeg and eg = —ée,., 


(11) 


which may be obtained by differentiating Equations (9) and applying the chain rule. 


Unit 28 Subsection 4.2 


The polar coordinates [r, 0] 
are defined in terms of the 

Cartesian coordinates (z, y) 
by the equations 


e-= reos#, 


g=rsind, 


Figure 8 


Equations (9) and (10) apply 
at all points of the plane 
other than the origin, where 
J,e, and eg are undefined. 


Figure 9 


These equations were derived 
in Unit 28 Subsection 4.2. 
Although attention was 
restricted there to circular 
motion, the derivations apply 
also in the general case. 
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We shall now apply Equations (11) to obtain expressions for the velocity v = r and 
acceleration a = r in terms of polar coordinates and unit vectors. As pointed out 
above, the position vector r is given by 


r=re,, (8) 
where r = |r|. Differentiating this equation and using the product rule, we have 

t=re,+reé,, 
which by the first of Equations (11) may be written as 

r= re, + réeg. ‘ers 
Differentiation of this expression for the velocity produces the relation 

i =e, +7re, + (764+ 1rb)ep + rbe, 


for the acceleration. After applying each of Equations (11) and collecting terms, we 
arrive at | 


i = (# — 6” )e, + (270+ rO)eg. 
Recognizing that 
ee 
276 +70 = —-—(r70 
— r dt ee, 
we can write the expression for the acceleration as 


‘ id ‘ 
eee eee tat Be 
r= (7 —r0“)e, + a (r“O)eg. 


The following box summarizes the results which we have obtained. 


(13) 


Position, velocity and acceleration in plane polar coordinates 


The plane polar unit vectors e,., eg are related to the Cartesian unit vectors 
i,j by the equations 


e, = cos#i-+ sin 6@j, €g = —sin#i+cos6j 


and i= cos0e, — sin@ eg, j=sin@e, + cos6eg. 


The derivatives e,,@g are given in terms of e,,eg by 


(11) 


In terms of r, 0, e, and eg, the position r, velocity r and acceleration ¢ of a 
particle moving in the (xz, y)-plane are given by the equations 


Y= re,. (8) 


r= re, + roeg, (12) 


Po ee aie, 


rit -) 


Of the following four exercises, the first three will give you practice in expressing a 
particle’s velocity and acceleration in plane polar coordinates. The fourth, Exercise 13, 
asks you to express the mechanical energy of a particle in polar coordinates. The result 
of this exercise will be applied later in the unit. 


Exercise 10 


A particle is constrained to move in a circle of fixed radius r = R whose centre is at the origin. 


Find the velocity r and acceleration r in terms of R, 6, 6, e, and eg. Give the velocity and 
acceleration also in terms of R, 0, 0, 0, i and j. 


Exercise 11 


A particle is constrained to move radially outwards from the origin, at some fixed angle 0 = 6p. 
Find the particle’s velocity r and acceleration ¢ in terms of 7, # and e,. Give the velocity and 
acceleration also in terms of 7, 7, 90, i and j. 


Note that 7 = d|r|/dt, the 
rate of change of the radial 
coordinate, is not in general 
the same as the speed |r| of 
the particle. 


These expressions for the 
velocity and the acceleration 
differ from those derived in 
Unit 28 Subsection 4.2 
because there we restricted 
ourselves to circular motion 
with r = constant. 


This is the particular type of 
planar motion which was 
considered in Unit 28 
Section 3 and Subsection 4.2. 
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Exercise 12 

A particle moves in the (z, y)-plane with radial coordinate r(t) = at” + bt and angle 6(t) = ct’. 
Find the velocity and acceleration in terms of t, e, and eg. 

Exercise 13 


It was shown in Subsection 1.2 that if a particle moves under the influence of a conservative 
force F = — grad U, then its total mechanical energy - = mr + U is constant. Show that if Recall that by r? we mean 
the motion takes place in the (z, y)-plane then FE may be written as rer = |r|’. 


EB =im(? +776?) +U. 
[Solutions on page 48 | 
Thus far we have derived kinematical formulas for the planar motion of a particle, 
using polar coordinates. Before proceeding to the subject of planetary orbits we shall 


briefly examine the mechanics involved in some simpler systems. The starting point for 
this examination is, as usual, Newton’s second law. 


In plane Cartesian coordinates, Newton’s second law mr = F takes the form 


m(Zi + ¥j) = F,i+ F,j, Here F, and Fy are 
; : Pas : respectively the z- and 
leading as in Unit 15 Section 2 to the two scalar equations y-components of the force F. 


me=—f, and my =F. 


Similarly, Newton’s second law mt = F may be expressed in plane polar coordinates, 
using Equation (13), as | 


; Lad : 
m (( — re, + = 5 (r78 00) = Fe, + Foes, Here F, and F% are 
r dt respectively the radial and 
giving the pair of scalar equations transverse components of the 
force F. 
m(#—r6é*) =F, and =a 8) mee (14) 
i 


It is these equations of motion which are sought in the following example and exercises, 
each of which represents a situation in which it is appropriate to use polar rather than 
Cartesian coordinates. 


Example 2 


A particle of mass m is launched at the bottom of a circular track of radius R (see 
Figure 11). As the particle moves along the track it is subjected to the downward force 
of gravity, a normal reaction force of magnitude F' perpendicular to the track, and a 
frictional force of magnitude f in the direction opposite to that of the particle’s motion. 


Set up the equations of motion for the particle in plane polar coordinates, assuming 
that the particle remains in contact with the track and moves in the sense of 
increasing 6. 


Solution 


The forces have been drawn on Figure 11, which also indicates our choice of Cartesian 
coordinate system. The reaction exerted by the track on the particle is —F'e,. The 
force of gravity is mgi, but i is given in terms of e, and eg by the first of 

Equations (10), so the gravitational force may be written as Figure 11 


mg(cos @e, — sin@eg@). 


The particle moves in the sense of increasing 0, that is, in the direction of eg, so the 
frictional force is —feg. Adding the three forces together gives the total force 


— Fe, + mg(cos@e, — sin eg) — fes = (-F' + mg cos @)e, + (—f — mgsin A)eg, 
so that the radial and transverse components of the total force are respectively 


F,=—F+mgcos@ and Fe = —f —mgsiné. 


14 


MST 204 30.1 


Since the particle remains in contact with the track, it is constrained to move in a circle 
with constant radial coordinate r = R, so that 7 = 7 = 0. The equations of motion 
(Equations (14)) are therefore 


—~ mR? = —F + mg cos 6, 
mRO = —f—mgsin@. O 


Exercise 14 


A puck of mass m moves without friction around the origin O on a horizontal table (see 
Figure 12). It is fixed to a light, inextensible string which is reeled in steadily, so that its 
distance from the origin is given by 


r(t)=r(0)—pt  (t< r(0)/p), 


where yp is a positive constant. The string remains taut throughout. 


(i) 


Show that the equations of motion for the particle which represents the puck are 


mré? =T and 5 (r70) = 0), 


where T’ is the tension in the string. 
Show that at time ¢ the tension is given by 
mr*(0)6?(0) _ mr*(0)6?(0) 
rs(t) (r(0) — pt)? 
Suppose that m = 0.1kg, 1» = 4 ms“, the initial distance is r(0) = 1m and the initial 


angular velocity is 6(0) = lrads~'. Suppose further that the string will break when it 
experiences a tension of 500 N. Find the distance r and the time t at which the string 
breaks. 


ues 


Exercise 15 


A puck of mass m moves frictionlessly around the origin on a horizontal table. It is connected 
to the origin by a perfect spring of natural length Jo and stiffness k, so that it is acted upon by 
a force 


(i) 

(ii) 
(iii) 
(iv) 


F = —k(r — lo)e,. 
Set up the equations of motion for the particle which represents the puck. 
Show that, however the particle moves, mr?6 is constant. 
Verify that a potential energy function for this force is U = sk(r — Io)”. 


Suppose that the puck has speed u when the spring has length 2/9. What is its speed 
when the spring has its natural length in the subsequent motion? 


[Solutions on page 48 | 


Summary of Section 1 


i. 


For particle motion in three or fewer dimensions, a force F acting on the particle is 
conservative if there exists a scalar field function U = U(r) = U(z, y, z) such that 


F= -— grad U. 


The function U is the potential energy of the particle. In such a case, the total 
mechanical energy of the particle, 


E=imr’?+U 


2 


(where r° = r-r = |r|*), is constant throughout the motion. 


In one dimension, the condition for a conservative force and the expression for the 
total mechanical energy become respectively 


d 
Pa" nd pine? oo 
dx 


where F' is the z-component of the force and U = U (a). 


With f = 0, these are also 
the equations of motion for a 
circular pendulum, where F' 
is now the tension in the 
string. You saw these 
equations derived via 
Cartesian coordinates in 
Unit 15 Section 5, and more 
directly in Unit 28 
Subsection 3.3. 


The upper bound on the time 
t is included because r(t) 
must always be non-negative. 
The puck would reach the 
origin at time t = r(0)/p if 
the situation remained as 
described until that time. 


Figure 12 


You will see this apparatus 
demonstrated and discussed 
further in the television 
programme for this unit 
(Section 4). 


MST204 30.2 15 


3. The unit vectors e,,e9 corresponding to the plane polar coordinates [r, 0] (where 
x =rcos@, y=rsin@) are related to the Cartesian unit vectors i,j by the 


equations 
e, = cos#i+sin6@j, €g = —sinfi+cos@j 
and i= cose, — sin@ eg, j=sinde, + cos eg. 


4. The derivatives €,,@g are given in terms of e,,eg9 by 
e,. = Jeg and eg a —e,. 


5. The position r, velocity r and acceleration r of a particle moving in the 
(x, y)-plane are given in terms of plane polar coordinates by the equations 


>= Te. 
r=re,+reg, 


ae ae LA 5: 
C= (7 3 r0" Je, =— ~ alr A)eo. 


2 Kepler’s laws and Newton’s law of 


gravitation 
The Earth is an ordinary-sized planet which orbits a middle-aged star of rather small 
size near the outer fringes of a typical spiral galaxy some 10° light-years across. This One light-year is the distance 
galaxy contains between 10'' and 10!” visible stars, together with much other matter, travelled by light in one year, 
and is one member of a rough association, or cluster, of some thousand galaxies of lee 
various shapes and sizes. 3 
Located in this myriad of astronomical objects is our solar system. It consists of a 
relatively massive Sun encircled by eight chief planets, some of which themselves Owing to its small size and 
possess one or more satellite moons. Any member of the solar system interacts via features of its orbit, Pluto is 


nowadays not considered to 


gravitational force with all of the others and, indeed, with all other heavenly bodies. 
be a planet. 


But distances between stars are measured in light-years—and the force of gravity 
between two bodies decreases rapidly in magnitude as their separation increases—so 
interactions with objects beyond the solar system can mostly be neglected. 


The gravitational influence of a body increases with mass, and within the solar system, 
there is a rough hierarchy of size. The Sun is by far the largest object, followed by the 
planets and then by their moons. The Earth’s mass is only about 3 x 107° that of the 
Sun, and its moon’s mass is about 10~ of its own. The largest planet, Jupiter, is 
roughly 300 times more massive than the Earth, but it is far away. In fact, the 
gravitational influence of a planet upon any other body is of little significance 
compared with the effect of the Sun, unless the body acted upon is relatively close to 
the planet. To a good first approximation, then, we can think of the planets as orbiting 
the Sun singly, and the moons as singly orbiting their planets. This is the model 
adopted by Newton, who considered the problem of only two bodies orbiting one 
another. (In a more advanced treatment, the effects of the other planets can be added 
as perturbations.) We shall assume in addition that one of the bodies is very much 
more massive than the other. 


An analysis of the motion of one relatively light body moving under the gravitational 
attraction of a second, relatively massive, body is the goal of Section 3. In this section 
we introduce the topic of planetary motion by describing (in Subsection 2.1) the 
features of their orbits which were observed by Kepler and stated in the form of three 
informative laws. These laws will be derived in Section 3 as a consequence of Newton’s 
universal law of gravitation, which we consider in Subsection 2.2. 
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2.1  Kepler’s laws described 


As they are worded, Kepler’s three laws refer specifically to the orbits of planets around 
the Sun, but they apply also to the orbit of any relatively light body about another 
much more massive body. The laws were established by Kepler from experimental 
observations of the motion of several of the planets. You will see later that they may be 
derived from Newton’s three laws together with his law of gravitation. It is the purpose 
of the current subsection to explain the meaning of Kepler’s laws. 


Kepler’s laws of planetary motion 
Law I Each planet moves in an ellipse, with the Sun at one focus. 


Law II The line joining a planet to the Sun sweeps out equal areas in equal 
times. 


Law III The square of the orbital period of a planet is proportional to the 
cube of the semi-major axis of its orbit. 


Kepler’s first law specifies the shape of a planetary orbit. Since it refers to an ellipse 
and to a focus, we start by considering the mathematics of ellipses. 


An ellipse is a symmetrical plane figure like that drawn in Figure 1. One of the 
implications of this first law, then, is that the orbit of a planet about the Sun lies in a 
plane. The usual definition of an ellipse is as a set of points (X,Y) which, with a 
suitable choice of origin and axes, satisfy the equation 


» ae 

a hae ae | (1) 
where a and 0b are positive constants with a > 6. The constants a and b are called 
respectively the semi-major axis and semi-minor azis of the ellipse. Their geometrical 
significance is shown in Figure 1. 


Exercise 1 


What can be said about an ellipse in the special case a = b? 


Exercise 2 


(i) Show that the area A of an ellipse can be expressed as 
A=2 [ 1 — (X/a)? dX. 


(ii) By using the substitution X = acosu, show that A = rab. 
[Solutions on page 49 | 


To understand Kepler’s first law, it remains to explain where a focus of an ellipse is 
located, and then to express the equation of the ellipse relative to this point. We shall 
express the equation with respect to the new set of axes Ozy shown in Figure 2 below. 
These axes are obtained from the set O’XY of Figure 1 by a translation of the axes to 
the right by a distance g, where 


q= Va? — b?. (2) 


(We could also shift axes to the left by the distance q, but our choice represents no loss 
of generality since this configuration is just the other rotated by an angle 7.) The 
points F’ and O of Figure 2, located a distance q to the left and to the right of the 
centre of the ellipse, are called the foci of the ellipse. 


Figure 1 


‘Foci’ is the plural of ‘focus’. 
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Figure 2 


Exercise 3 
Where are the foci of a circle? 


(‘Solution on page 49 | 


Since the translation is along the coincident «- and X-axes, the new coordinates (2, y) 
are related to the original ones (X,Y) by 


A= 2+, =e 
Then, from Equation (1), the equation of the ellipse takes the form 
Cee 
a ee 
We shall now proceed to re-express this equation in plane polar coordinates. After 
multiplying the equation by b? and rearranging, we obtain 


ae 


b b? 
y = b°— J(a +9)" = 0° - G (a + 2qr +9"). 


Adding x? to both sides and collecting terms gives 
b? 2b*q b? 
x + y? — ( mee a”) pre eS eae (: = 7) x, 


The left-hand side is r?, the square of the radial distance, so that, on using 
Equation (2) to express qg in terms of a,b, we have 


b4 bx b? b? 
2 2 
re = —~ —2—4/1-—-—- + (1-4) 2’. 3 
a? a a? ( a (3) 
This rather cumbersome expression is made neater by employing two new quantities 
defined in terms of a and b. These quantities are the eccentricity This usage of the symbol e 
should not be confused with 
eee b? 4 its use as the base of natural 
ae (4) logarithms, 2.718 28... 
and the semi-latus rectum 
b 
[=—. (5) 
a 


Note that, since b < a, we have 0 < e < 1. This is a condition on the eccentricity for 
ellipses. 


Exercise 4 

Find e when 

(i)a=b (acircle); (ii) b = Sa; (iii) b = $a. 

Comment on the relationship between the eccentricity and shape of an ellipse. 


[Solution on page 49 | 
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From Equations (3)—(5), we now have 
ee ge e247? 
or ee es a 
Then taking the square root of both sides gives 
Pee tl ee), (6) 


We must now decide on the choice of sign for the right-hand side of Equation (6), to be 
consistent with the fact that r is never negative. We shall show that, for all points z on 
the ellipse, the quantity | — ez is positive. From Figure 2, the largest value of x on the 
ellipse is « = a —q. Thus we have, for all points on the curve, 


!—ex >1—e(a—4q). 


But we also have gq = Va? — b?, | = b*/a and e = q/a, so the right-hand side of this 
inequality is 


| 2 Eee: 
|—e(a—q)= ~~ Y*—* (q_ a? — b) 
a a 
2 2 2 
se a2 gies 
a 
=a-— Va’ — Bb, 


which is never negative. Hence / — ez itself is never negative, when x assumes values for 
points on the curve. Thus the positive sign is appropriate in Equation (6), and the 
equation for an ellipse with respect to the axes Oxy of Figure 2 is 


l 
r=1—exr -or ~ tte 
3 ir 
Since x = rcos@ (see Figure 3) we have the convenient form 
l 
-=1+ecos6 (7) 
. 


for the equation of an ellipse in polar coordinates with origin at a focus of the ellipse. 


Kepler’s first law therefore says that if the origin of coordinates is located at the Sun, 
then the orbit of each planet lies in a plane containing that origin, and is described by 
Equation (7) once the axes have been appropriately oriented. With a different 


orientation for the axes (but the same origin), Equation (7) becomes 
l 
— =1+ecos(6 — 6), (8) 
i 


where 6 is the clockwise angle through which the axes have been rotated from their 


position in Figure 3. The major axis of the ellipse then lies along the direction specified 
by 6 = Oo. 


Kepler’s second law says that the line joining a planet to the Sun sweeps out equal 
areas in equal times. This law can be used to specify the rate at which a planet moves 
in its particular elliptical orbit. 


Suppose, without loss of generality, that the planet moves in an anticlockwise sense. In 
Figure 4 the two shaded regions represent the areas swept out by the planet between 
times ¢; and t2, and between times ts and ty. Assuming that the two time intervals are 
of equal duration, Kepler’s second law says that these areas are equal. If A(t) is the 
total area swept out since some initial time, t = 0 say, this law can be expressed as 


A(t4) = A(t3) = A(t2) oor A(t,) whenever t4 at t3 = to a ty. 


Putting xz = 0 into this 
equation gives r = +l, 
showing that the semi-latus 
rectum / is the distance from 
the origin of the ellipse’s 
intercepts on the y-axis. This 
is illustrated in Figure 3 
below. 


Figure 3 
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Figure 4 


A more concise way of stating Kepler’s second law is that the rate at which area is 
swept out is constant, that is, 


A is constant. 


It is possible also to state the law in a form which relates the radial coordinate r to the. 
rate of change @ of the angular coordinate. Figure 5 shows a close-up of part of the 
planetary orbit, where we consider the area 6A swept out during a small time interval 
of duration 6t. Over this interval the planet moves from point P at time t to point P’ 


at time t + 6t, traversing an angle 66 and changing its radius from r to r + 6r. If ét is poe 
sufficiently small then the area 6A will lie between the areas of two circular sectors with 
angle 60, one having radius r and the other radius r + ér. If r(t) is increasing in the If r(t) is decreasing then the 
part of the curve under consideration then this gives expressions for the upper and 
lower bounds in this 
7760 <é6A< $(r br)? 60. inequality are reversed, 


oe 5 coe leading to the same outcome. 
On dividing through by é6t, and then taking the limit as é6¢ tends to zero, we obtain 


A = ir. (9) 


Since A is constant according to Kepler’s second law, this law may also be interpreted 
as saying that 


r?@ is constant. 


Equation (9) can also be used to express an area swept out by the planet in terms of r 
and 9. Thus the shaded area in Figure 4 covered between times ft, and f2 is 


A(tz) — A(t) = [4]? = / 


Exercise 5 


A body moves in a circular orbit around the origin, with radius R and fixed speed V. Show 
that A is constant, and find its value. 


Exercise 6 


Compare the average speeds of the orbiting planet during the two equal intervals t; < t < te 
and t3 <t < t4 in Figure 4. Your answer should consist of a few sentences only. 


[Solutions on page 49 | 


Kepler’s third law gives a relation between the period of the planetary orbit and the 
lateral extent (semi-major axis) of the orbit. If T is the period then, by definition, the 
planet will take this time to make one full circuit of the ellipse. If the particular orbit 
has semi-major axis a, then Kepler’s third law says that T? is proportional to a’, 

that is, 


T? = ka’, 


where k is a constant. This law, then, makes a statement connecting the overall rate of 
motion of the planet to a geometrical property of its orbit. We shall prove this result, 
together with Kepler’s other two laws, in Section 3. 
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Exercise 7 


The semi-major axis of the planet Mars is 1.524 times as large as the Earth’s semi-major axis. 
The Earth’s period for its orbit around the Sun is 365.256 days. Find the period (in Earth 
days) for the orbit of Mars around the Sun. 


[Solution on page 49 | 


2.2 Gravitational force and potential energy 


I began to think of gravity extending to ye orb of the Moon, and... from 
Kepler’s Rule ... I deduced that the forces which keep the Planets in their 
Orbs must [vary] reciprocally as the square of their distances from the centres 
about which they revolve: and thereby compared the force requisite to keep 

_ the Moon in her Orb with the force of gravity at the surface of the Earth, and 
found them to answer pretty nearly. All this was in the two plague years of 
1665 and 1666, for in those days I was in the prime of my age for invention, 
and minded Mathematicks and Philosophy more than at any time since. 


Isaac Newton’s account of his discovery of the universal law of gravitation. 


This subsection starts with the statement of Newton’s universal law of gravitation, then You met Newton’s law of 
derives the potential energy function corresponding to the (conservative) gravitational gravitation in the television 
force exerted by a particular type of extended object. Finally, we compare the theory of Programme for Unit 28. 
gravitation introduced here with the simplified version which has been applied in 

previous units for motion near the surface of the Earth. 


Newton’s universal law of gravitation describes the nature of the force between pairs of 
particles. It gives the attractive force between them as a function of the masses of the 
particles and of their separation. It is universal in the sense that it is hypothesized to 
hold between the particles comprising all bodies in the universe, and not just between 
those comprising the Earth and an apple, or the Earth and the Moon. In words, the 
law can be stated as follows. 


Newton’s universal law of gravitation 


The force of gravity between two particles is attractive, directed along the line 


between them, proportional to the product of their masses, and inversely 
proportional to the square of their separation. 


Let us write this mathematically, using the vector notation which was unavailable to 
Newton. We write the position vectors of the two particles as r;,r; and their masses as 


m;,m;. The gravitational force on m; due to the presence of m,; is denoted by ry. As a shorthand device, we 
From the verbal statement of Newton’s law above, we know the following. shall refer to the particles by 
using their masses m;,m, as 
(i) Fy; is directed from m; to mj, hence its direction is defined by the unit vector labels. 
rT . = 
j 


ae (see Figure 6). 
g—T% 


1 


(ii) The magnitude of F,; is proportional to mjim,;, and to | Uh 
Be £y 


Putting these elements together gives the force of gravity on m; exerted by m; as 
Gmim,; 


Fi; = ———4(r; — ri). | 
J oe r;) (10) 


See ; wre O 
This is Newton’s universal law of gravitation in vector form. Here G is the 


constant of proportionality, called the universal constant of gravitation, or gravitational Figure 6 
constant. It is one of the least precisely measured physical constants. Its value, in SI 
units, is 


G= 647 x10". (11) 
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Exercise 8 


What are the SI units of G? 


Exercise 9 
Yr; ae! i 


Show that is a unit vector. 


lr; —r i| 
Exercise 10 


Gravitation is an inter-particle force. By Newton’s third law (Unit 17 Section 2), the force 
exerted upon m; by m; should be equal in magnitude, but opposite in direction, to the force 
exerted upon m; by m;. Show that the universal law of gravitation is consistent with this, so 
that ei; + PA == 3. 


[Solutions on page 49 | 


Most of the rest of this subsection is taken up with consideration of the gravitational 
attraction caused by an extended object with spherically symmetric mass distribution, 
and calculation of the corresponding potential energy function. 


Using Equations (10) and (11), it is possible to calculate the force of attraction between 
any two objects, provided that they are sufficiently small to be modelled by particles. 
In many cases, however, it is far from the truth to claim that both objects are ‘small’ 
relative to the situation being considered. For example, Figure 7(a) shows a satellite of 
mass m, in orbit around the Earth at a height approximately equal to the Earth’s own 
radius. At first sight it seems unreasonable here to model the Earth by a particle, 
although it appears sensible (for most purposes) to use a particle model for the satellite. 


(a) (b) 


Figure 7 


In fact, you will see in this subsection that if the Earth’s mass is assumed to be 
distributed in a spherically symmetric manner (to be defined below), then its attraction 
on the satellite is the same as if all its mass were concentrated in a particle at its centre 
(see Figure 7(b)). 


To approach this problem, consider the gravitational force on some particle of mass m 
outside an arbitrary extended object. We want to calculate the net force acting on the 
particle with respect to some coordinate system Oxvyz, as in Figure 8. To do this, we 
regard the extended body as being composed of many small volume elements. Then the 
net force on the exterior particle is obtained by adding together all the (vector) 
contributions from these volume elements. Since the elements are small, they can be 
treated individually as particles. 


If you are short of time, you 
might like to read as far as 
the paragraph below Figure 7 
and then move directly to the 
summary box on page 25: 
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Figure 8 


The force exerted upon the particle of mass m by a typical volume element, which has 
mass m; and position vector r;, is obtained by applying Equation (10). This force is 
pea 3 


Gmm; ae 


Ir 
where r is the position vector of the mass m. Adding together the contributions from all 


the volume elements in the extended object gives the net gravitational force on m as This is the gravitational force 
— on the particle of mass m 
mi(ri — Yr) due to the system of N 
F(r)=Gm)_ ame (12) particles. Although this is 
= es only an approximation to the 


. gravitational effect of the 
By construction, the sum in Equation (12) supposes that the extended ob ject consists Pleated objeck ike 


of a large number, N, of small elements m;. We shall ultimately consider (for a approximation becomes exact 
particular case) the limit of this sum as N becomes large and each m; tends to zero, in the limit as N is increased 
which leads to an expression for F(r) as an integral over the volume of the extended indefinitely. 

object. This could be done by taking the limit directly from Equation (12), but it is 

simpler to rewrite it first in terms of gravitational potential energy. 


You saw in Subsection 1.2 that a force F(r) is conservative, with potential energy 
function U(r), provided that 


F = — grad U. Equation (3) of Section 1 


ae (page 9) 
We shall now show that the net force given by Equation (12) can be expressed in this 
way. The main step in this demonstration is to establish that 


| ae 


= grad (13) 


lr — r;|° Ir —r;| 


In order to prove Equation (13), it suffices to use Cartesian coordinates, in terms of 
which we have 


for = (x Be ae a oe yi)” ie A ee is 
The definition of the gradient gives 
2 1 23 1 O 1 
r—r| ‘Oclr—r] Sdyfr—n] | Oz [r—e 
The i-component on the right-hand side of this equation is 


2 fee 
Ox |r—r;| 


= i(—3) (( — 4)? + (y—:)? + (2 — 4)?)” x (ae — 2) 


grad (14) 


Oe ae as 


ee 
Ir —r;|° 


The calculations for the j- and k-components of Equation ( 14) are similar. Using the 
relations 


r=zi+yj+zk and rj = xj41+ yj t+ zk, 
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we obtain 


grad = ie oly tte a) os 


i eae 


lr —r;| 


a lr — r;|?" 
which proves the validity of Equation (13). 
Combining Equations (12) and (13) then gives the gravitational force F in terms of a 
potential energy function, as 


N 
1 
Fir} = Gm S_m; grad : 


N 
mj 
= grad | Gm 
Eocene a) 
or P(e) = era (- om). 


This equation expresses the net gravitational force exerted by a system of N particles, 
with masses m, and corresponding position vectors r; (¢ = 1,2,...,N), upon an 
external particle of mass m and position vector r. It has the form 


F(r) = — grad U(r), 


where the potential energy function U is given by 


U(r) = —Gmy° —. (15) 


Exercise 11 


A particle of mass M is fixed at the origin. Find the corresponding gravitational potential 
energy function U for a second particle which has mass m and position vector r. Find also the 
gravitational force exerted upon the second particle by the first. 


[Solution on page 49 | 


We seek now to progress from the gravitational potential energy of a particle due to a 
system of particles to the corresponding potential energy for a continuous distribution 
of mass, as in Figure 8. To achieve this, we characterize the distribution of mass by its 
mass density function, p(r). This function is defined so that any small volume element The mass density function p 


having volume 6V; and centred at position r; contains a mass m;, where was introduced in Unit 27 
Section 2, where it was used 
m; ~ p(ri)6V;. to calculate the overall mass 


: ae : < of extended bodies. 
Inserting this expression into Equation (15) gives, for the gravitational potential energy 


U of a particle due to a body with continuous mass distribution, 


sa yo seat a 


lr —r;| 


This approximation to U becomes exact in the limit as the size of the volume elements 
6V; tends towards zero and their number, N, increases indefinitely. In this limit, we 
obtain 


ae 


big oa , AE) ay | (16) 


The right-hand side here is a volume integral with respect to the variable r’, extending 
over the entire body B. 


Equation (16) gives quite generally the gravitational potential energy of a particle of 
mass m located outside any continuous distribution of mass, specified by the mass 
density function p. We now specialize to the case of a mass m lying outside a 
spherically symmetric distribution of mass. 


By definition, a mass distribution is spherically symmetric if the corresponding mass 
density function p(r) depends only upon the radial distance, r = ,\/x? + y? + z?. We 
express this spherical dependence of the mass density function by writing it as 


p(r) = p(lr|) = e(r). 
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Exercise 12 


Which of the following mass density functions p correspond to spherically symmetric mass 
distributions? (A, B and R are positive constants.) 


(ay “p(r) = A ie ty? + z) esa Cas 
(b) p(r) = A (a? +y?) aes (r > 0) 
ic) ey =A ie +y? + ey Pe a (¢ >) 
(d) p(r) = fg : : R) es 
a ee 

A <0 S9-= FR) 


(f) p(r) = ‘6 (x > RY. 


[Solution on page 49 | 


The integral on the right-hand side of Equation (16) was considered in Unit 27 for the 
case of a spherically symmetric mass distribution. It was proved there that when 


p(t’) = p(|r’|) then 


/ 
Ar) ay il (17) 
p\r’-r| r 
where 
M= | p(r’) dV (18) 
B 


is the total mass of the sphere B. By Equations (16), (17) and (18), the gravitational 
potential energy function for a spherically symmetric mass distribution is 
GmM 


U(r) = - as 


(19) 


Exercise 13 
Show that the gravitational force exerted on a particle of mass m due to a spherically 
symmetric mass distribution of mass M with its centre at the origin is 


GmM 


p2 


F(r) = — e,, 


where e, = r/|r| is the unit vector in the outward radial direction. 


[Solution on page 49 | 


Equation (19), together with the results of Exercises 11 and 13, establishes the truth of 
the statement following Figure 7: the gravitational attraction of a body with spherically 
symmetric mass distribution on an exterior particle is the same as if all of the body’s 
mass were concentrated at its centre. In proving this result we have derived expressions 
for the gravitational force and potential energy which will be applied in Section 3 to the 
task of specifying planetary orbits. These expressions are given once again below. 


Unit 27 Subsection 4.3 


If the body is of finite extent 
and has a spherically 

symmetric mass distribution, 
then it is necessarily a sphere. 
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The gravitational effect of a spherically symmetric body 


A body has a spherically symmetric mass distribution if (with the origin of 
coordinates at the centre of the body) its mass density function p(r) depends 


only on the radial distance |r| = r = ,/x? + y? + z?. The gravitational force 
exerted by such a body B on an exterior particle with mass m and position 


vector r is given by | It is said in some quarters 
Con Mt that it took idee twenty 
ae e,., years to prove this result, 
- ; causing him to delay 
where e, = r/r is the unit vector in the outward radial direction, publishing Principia. But 
this version of events has not 
GmM been established to the 
Sas satisfaction of all science 
historians. 


F(r) = — grad U(r) = 


is the gravitational potential energy of the particle, and 


M= | o(e')av 


is the total mass of the body. 


Exercise 14 


Two spherically symmetric bodies have masses M and m. If their centres are separated at 
some instant by a distance R, what are then the magnitude and direction of the gravitational 
forces exerted by each body upon the other? {Hint: By a result of Unit 17 Section 2, the 
motion of the centre of mass of a body is the same as that of a particle of the same total mass 
which experiences all of the external forces applied to particles of the body.| 


[Solution on page 50 | 


To conclude this section, we examine how Newton’s universal law of gravitation matches 
up with the simpler model used previously for motion near the surface of the Earth. 
When modelling motion close to the Earth it is customary to assume, as we have done 
in earlier units, that the force on a particle of mass m is directed vertically downwards 
and has the constant magnitude mg, where g = 9.81ms ”. More generally, under the 
assumption that the Earth is a spherically symmetric distribution of mass, we have just 
shown that at any point on or above the Earth’s surface the gravitational force is 
mMG 


F(r) = —- a ae 


where M is the mass of the Earth, e, = r/r is the unit vector in the outward radial 
direction and r = |r| is the distance from the centre. 


Clearly the two descriptions agree as to the direction of the gravitational force, since 
—e,., being directed towards the centre of the Earth, is always ‘vertically downwards’ 
from a local point of view. It remains then to compare the two expressions for the 
magnitude. 


Exercise 15 
(i) Show that g = MG/R?, where R is the radius of the Earth. 


(ii) Given that g = 9.81ms~*, G = 6.67 x 10° m?kg"'s-? and R= 6.4 x 10° m, estimate 
the mass M of the Earth. 


[Solution on page 50 | 


It can be shown that the model of a constant gravitational force is quite a reasonable 
one for heights of a few kilometres above the Earth’s surface. Consider two points P, 
and FP, lying on the same radial line, with P,; on the surface and P 2 at a distance h 
above it (see Figure 9). We wish to compare the force of gravity at points P,; and P, on 
a particle of mass m. 


At both points the force on m is ‘downward’, along the direction of —e,. The 
magnitude of the force at P; on the Earth’s surface (at radius R) is 
mMG 

Re 


Figure 9 
Fy= 
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At the point P2, which is a height h above P,, the force has magnitude 


oe mMG 
The ratio of these two magnitudes is 
Fy R? 1 


————_—___ 


Fee a/R 


Provided that h is small compared to R, the force of gravity at point P> is nearly equal 
to that at P,. 


Exercise 16 


Given that the acceleration due to gravity is 9.81ms~? at the Earth’s surface, and that the 
radius of the Earth is 6.4 x 10° m, estimate the acceleration due to. gravity at a height of 10*m 
above the Earth’s surface. 


[Solution on page 50 | 


Summary of Section 2 


1. Kepler’s laws of planetary motion are as follows. 
(I) Each planet moves in an ellipse, with the Sun at one focus. 
(II) The line joining a planet to the Sun Sweeps out equal areas in equal times. 


(III) The square of the orbital period of a planet is proportional to the cube of the 
semi-major axis of its orbit. 


2. (i) The ellipse shown in Figure 10 has semi-major axis a and semi-minor azis b. 
With respect to the axes O’XY, its equation is 


xX? y2 
—— =], 


= 


where a > b> 0. 


Figure 10 


(ii) The area of the ellipse is ab. 


(iii) The foci of the ellipse are located a distance gq to either side of the centre O' : 
where 


q= V a? — B?. 


(iv) With respect to axes Oxy with origin at the right-hand focus, and in polar 
coordinates, the equation of the ellipse is 


l 
—-=1+ecos8@, 
T: 


MST204 30.3 27 


where 
b 


a 


l 


is the semi-latus rectum, and 


is the eccentricity, where 0 < e < 1. 


3. Newton’s universal law of gravitation states that the gravitational force on a 
particle with mass m; and position vector r; due to another particle with mass m; 
and position vector r; is 


(r; a t;), 


where G is the gravitational constant. 


Fr. - Gmim, 
so 


4. A body has a spherically symmetric mass distribution if (with the origin of 
coordinates at the centre of the body) its mass density function p(r) depends only 
on the radial distance |r| = r = \/x? + y? + z?. The gravitational force exerted by 


such a body B on an exterior particle with mass m and position vector r is given 


by 
GmM 


r2 


F(r) = — grad U(r) = — 


er, 


where e, = r/r is the unit vector in the outward radial direction, 


Z GmM 


Ce = 


fs 


is the gravitational potential energy of the particle, and 


M= i. alae haV. 


is the total mass of the body. 


> Orbis 


In Subsection 2.2 Newton’s universal law of gravitation was introduced, leading to a 
derivation of the corresponding gravitational potential energy function for a spherically 
symmetric mass distribution. In the current section we shall use this result in 
developing a mathematical model to describe the orbit of a body under the 
gravitational attraction of a much larger body. To a reasonably good approximation, 
this model applies to the orbits of a planet around the Sun, of a moon around a planet, 
of a man-made satellite around a planet or around the Sun, and of comets around the 
Sun. The model will enable us to derive Kepler’s three laws of planetary motion, which 
were stated and explained in Subsection 2.1. 


The modelling assumptions which suffice to yield these laws are as follows. 


Modelling assumptions for gravitational orbits 


(i) The world-view of Newton is valid, including his three laws of mechanics 
and his universal law of gravitation. 


All forces other than the gravitational attraction between the two bodies 


may be neglected. (In particular, the gravitational effect of other 
celestial objects is ignored.) 


One body is much more massive that the other. 


The more massive body has a spherically symmetric distribution of mass. 


Kepler’s laws are also the 
first item in the Summary of 
Section 2 on the previous 


page. 
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Of these assumptions, (i) is the hardest to improve upon and (iii) is the easiest, 

while (ii) and (iv) are intermediate. It was only with Einstein’s general theory of 
relativity in the early twentieth century that small corrections to Newton’s world-view 
(Assumption (i)) were introduced. This was done in order to explain small departures 
from the Newtonian prediction that had been observed in the orbit of the planet 
Mercury. We could (but will not!) improve upon Assumption (iii) using the 
mathematics developed in the course. This improvement is required, for example, to 
describe the orbits of two stars of comparable mass under their mutual gravitational 
attraction. 2 


3.1 Central forces 


By Assumption (iii) above, we suppose that the more massive body (of mass M) does 
not move at all under the gravitational action of the lighter body (of mass m). Its 


centre of mass coincides with its geometric centre, by Assumption (iv). For convenience, 


we choose the origin of coordinates to be at this centre, as shown in Figure 1. 


Figure 1 


It was shown in Unit 17 that the centre of mass of a system of particles acted upon by 
external forces moves as if all the mass of the system were concentrated there. We may 
therefore treat the motion of the centre of mass of the lighter body as if it were a 
particle. Furthermore, since the heavier body is assumed to be spherically symmetric 
(in gravitational terms), we can apply the result of Subsection 2.2 which equates the 
gravitational effect of such a body with that of a particle of mass M at its centre. The 
resulting two-particle view is shown in Figure 2, where r is the position vector of the 
lighter particle. - : 


There is a simplifying feature of this model which can be shown at the outset: the 
motion of m occurs in a plane. The proof of this statement rests on the fact that the 
gravitational force exerted on m is a central force. In general, this term describes a 
force which is directed from the particle m towards or away from the source of the 
force. If the source is at the origin, this means that a central force has the direction of 
—r or of r. The definition is as follows. 


A central force F, acting from the origin on a particle with position vector 
r, has the form 


Fir) = oir, (1) 


where g(r) = g(x,y, z) is a scalar function of position. 


Exercise 1 


Which of the following are central forces (where K is a constant)? 


(a) F(t) =Kr (a) P(r) = +33 
(bh) Fig “ ie Fin = = + j4-k) 
<— S aor e tao= ~ (ai +a} + 2k) 


Unit 17 Subsection 2.3 


A similar argument to this 
was employed in the solution 
to Exercise 14 in Section 2. 


Equation (19) of Section 2. 


As earlier in the unit, we use 
the symbol m as a convenient 
label for the particle of mass 


re 


Figure 2 
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Exercise 2 
What is the scalar function g(r) for the gravitational system of Figure 2? 


[Solutions on page 50 | 


The angular momentum L of a particle with mass m and position vector r relative to 
the origin is, as defined in Unit 29, Unit 29 Subsection 1.2 


Lax on. (2) In Unit 29 we used the 
symbol l for the angular 


This definition is not confined to two-dimensional motion, though the main application omentum of a single 


in Unit 28 was to motion in a circle. However, it is a consequence of the results of particle and L for the angular 
Exercise 3 below that if the particle is acted upon by a central force then its motion momentum of a system of 
does lie in a plane. particles. Throughout this 


unit we refer only to the 
angular momentum of a 
single particle, and for visual 
Show that: convenience this is 
represented by the symbol L. 


Exercise 3 


(i) when the anguJjar momentum L is non-zero, the position vector r is always perpendicular 
to L; ; 


(ii) when the force on the particle of mass m is a central force then L = 0, so that L is a 
constant vector; 


(iii) when L = 0, the velocity r and position vector r are parallel or anti-parallel if neither is 
zero. 


[Solution on page 50 | 


From parts (i) and (ii) of Exercise 3, we can say that if L 4 0 then the position vector r 
of the particle of mass m is always perpendicular to the fixed direction defined by the 
constant angular momentum vector L. In the general case then, when L is non-zero, 
the particle acted upon by a central force moves in the plane perpendicular to the fixed 
direction of L, as shown in Figure 3. Also (from part (iii)), in the event that L is zero, 
the velocity lies along a straight line from the origin, and the particle will continue to 
move along this straight line. 


Figure 3 


The constancy of the angulay momentum L for a particle in gravitational orbit arises 
from the fact that gravity is a central force, and not from any other detail of its form. 
This result applies also, therefore, to other central forces. 


Motion under a central force 


A particle with mass m and position vector r, which moves solely under the 
action of a central force, has constant angular momentum 


L=f * mr. (2) 


Its motion is confined to the plane perpendicular to the direction of L when L 
is non-zero and to a straight line if L = 0. 
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As the motion is planar, we can now choose to orient the coordinate axes so that the 
motion lies in the (x, y)-plane, as indicated in Figure 4. This will permit us to apply 


the apparatus of plane polar coordinates which was developed in Subsection 1.3. It has yet to be established 
that the planar motion of a 
particle in gravitational orbit 
is as simple as that depicted 
here. For example, some 
other central forces lead to 
spiralling motions. 


Figure 4 


Since the angular momentum L is a constant vector, both of its direction and 
magnitude are constant. The constancy of its direction sufficed above to show that the 
orbit always lies in a plane. This is actually part of the content of Kepler’s first law, 
which asserts that the orbits are ellipses and therefore plane figures. 


We have not yet used the additional information that the magnitude of L is constant, 
and this is actually sufficient to prove Kepler’s second law. As you may recall from 
Subsection 2.1, this law relates to the rate at which area A is swept out by the line 
joining a planet to the Sun (which corresponds here to the line of the position vector r). 
We showed that A was related to the radial distance r and the angular velocity - by 
the equation 


A= 176. Equation (9) of Section 2 


Kepler’s second law states that A is constant, which from above is equivalent to the 
condition 


© (736) =0. (3) 


In the exercise below you are asked to show that Kepler’s second law is a consequence 
of the constancy of angular momentum. 
Exercise 4 


(i) By expressing r and fr in terms of the plane polar unit vectors e, and eg (Equations (8) 
and (12) of Section 1), show that the angular momentum of a particle moving in the 
(x, y)-plane is L = L.k, where 


Ly = mr°6. 


(ii) Deduce that if L, = 0 then Equation (3) is satisfied. 
[Solution on page 50 | 


In summary, thus far, we can say the following. 


Kepler’s second law, and the fact that each planet moves in a plane, are 
consequences only of the fact that the gravitational force between the planet 


and the Sun is a central force; taken together, they amount to a statement of 
the constancy of angular momentum. 
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3.2 Isotropic central forces and gravity 


With the force of gravity in mind as the ultimate target, we now specialize from central 
forces to the more restricted case of isotropic central forces. These are central forces for 
which the magnitude depends only on the radial distance from the origin. 


An isotropic central force has the form 


F(r) = f(r)e,, (4) 


where r = \/x* + y? + 2? is the distance from the origin, and e, = r/r. 


Exercise 5 
(i) Which of the forces listed in Exercise 1 (page 28) are isotropic central forces? 
(ii) Show that the gravitational force is isotropic, and find the function f(r) for this case. 


[Solution on page 50 | 


Choosing our axes as in Figure 4, so that the orbit lies in the (x, y)-plane, we have 


z=0Oandr= \/x?+y?. Newton’s second law, for the motion of a particle of mass m 
acted upon by an isotropic central force, is 


mr = {(rje,. 


Since the right-hand side ofthis equation is directed along the unit vector e,., it makes 
sense to express the acceleration r in terms of plane polar coordinates, using 

Equation (13) of Section 1. The particle’s equation of motion in the plane can then be 
written as 


oe m d a he 
m(F — ré“)e, + E a O)eg = f(rje,, 


from which we have 


m(# — 76?) = f(r) (5) 
™m d 2 Be cae 
and = 0) = 


The second equation is equivalent to Equation (3), whose validity was established for 
any central force in Exercise 4. The consequent constancy of the angular momentum 
magnitude, 


L. = mr76, (6) 


can be exploited to eliminate the variable 6 in Equation (5). Thus on putting 
9 = L,/mr? into Equation (5), we obtain 


or mr = f(r) + —. (7) 


Equations (6) and (7) contain, between them, the possibility of a complete solution to 
the problem of finding the orbit of a particle of mass m under the action of any isotropic 
central force, given by Equation (4). The procedure to do this would be as follows. 


(i) The initial conditions are the values of [r, 6] and [7,6] at some initial time, say 
t= 0. 


(ii) Since L, is a constant, these initial conditions determine the value of L, for all 
subsequent times, from Equation (6). 


(ili) Equation (7) may be solved, either analytically if possible or numerically, to obtain 
an expression or values for r(t), where t > 0. 


(iv) This solution may be substituted into Equation (6) to obtain an expresion for 0(t), 
from which 6(t) can be found by direct integration. 


A central force was defined 
by F(r) = g(r)r, which is 
equivalent to 

F(r) = g(r)re,. 
This force is also isotropic if 
g(r) = g(r). 
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The approach we shall eventually adopt, however, is somewhat different. Instead of 


seeking r(t) and @(t) individually, which are parametric equations for the orbit, we look 


for a direct relationship r = r(@). This may describe an orbit which is bound or 
unbound. For a bound orbit (see curve (a) in Figure 5), the distance of the particle 
from the origin O is bounded for all time. For an unbound orbit (see curve (b) in 
Figure 5), the particle approaches from infinitely far away, is deflected by the force 
centred at O and departs again to infinite distances. Corresponding to the bound orbit 
shown in Figure 6(a), the function r(@) has the graph of Figure 6(b). 


Figure 6 


Thus far, the results obtained relate to motion under any isotropic central force, and 
not just to gravity. At this point we shall finally specialize to the case of gravitation, 
though the results may also be applied to any other force which is described by a 
similar inverse square law. 


Specializing to gravity means using the fact that for this particular isotropic central 
force we have (from Exercise 5(ii)) 


GmM 
f(r) ae Sal r2 ’ 
so that Equation (7) for the particle motion becomes 
i? Gav 
me i g 
ae r2 (8) 


This equation looks somewhat formidable, but it can be tamed. In Subsection 3.3 you 
will see how a relationship r = r(@) can be obtained from Equation (8) to describe 
gravitational orbits. 


For the remainder of this subsection we shall investigate what can be deduced about 
the motion without solving Equation (8) directly. This indirect approach depends on 
the constancy of both the angular momentum L, and the total mechanical energy, E. 


It follows from the result of Exercise 13 of Section 1 that, for a particle moving in two 
dimensions and acted upon by the (conservative) gravitational force, the constant total 
mechanical energy E may be written as 


sm(r* + 7767) + U(r), | (9) 


where U(r) = —-GmM/r is the gravitational potential energy function. From 
Equation (6), the constant angular momentum is given by L, = mr?6. After using this 
equation to replace 0 by L,/(mr?), Equation (9) becomes 


2 
E = 5m + kmr? ( = | —- 
mr r 
L? GmM 
S es (sc a : ) 


Hence, by using the constancy of the angular momentum, we have expressed the total 
mechanical energy in terms of 7 and r only. 


Figure 5 


The electrostatic force is one 
example of an isotropic 
central force. You considered 
a two-dimensional example of 
such a force in Exercise 15 of 
Section 1. 


Equation (19) of Section 2 


We write U(r) rather than 
U(r) here to emphasize that 
the gravitational potential 
energy U depends only upon 
the radial distance r. 


In fact, any isotropic central 
force is conservative, and has 
a constant total mechanical 
energy F given by 

Equation (9). You will see 
this shown in Section 4. 
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As a result, Equation (10) has the appearance of an energy equation for a particle As a check, if Equation (10) 


moving in one dimension. The first term, $m”, resembles the one-dimensional kinetic _ is differentiated with respect 
energy. Since the bracketed term in Equation (10) looks like a one-dimensional to time (with F = 0), then 
the outcome is Equation (8). 


ena energy function, we shall call it the effective potential energy and denote it by Pad-auotlee ey, 
Ue"(r). This gives Equation (10) is a first 
.2 f integral of Equation (8). 

E = 5mr* +U(r), (8) 

where the effective potential energy is 
Se GmM 
U(r) = —*, - 7 (11) 
Dir? 


In summary, then, the radzal motion of the particle is as if it were moving in one 

dimension, with r > 0, under the influence of a one-dimensional conservative force 
having potential energy function U*"(r). From this insight you might expect that 
Equation (8) for the radial motion could be written as 


d 
mi = ——U*(r), 
dr 2 
by analogy with the definition of a one-dimensional potential energy function. This is Equation (1) of Section 1 
indeed the case, as can be verified by differentiating Equation (11). 


Obviously, the particle whose orbit we seek is not really moving in a single dimension, 
since as r varies with time the angle 6 is changing too. However, we have managed in 
Equation (10) to express the total mechanical energy of the orbiting particle in the 
same format as that for a one-dimensional motion. This means that we can deduce 
information about the radial coordinate of the orbit, using the effective potential 
energy function, with the approach that was adopted in Subsection 1.1 for an actual 
one-dimensional motion. 


You saw in the analysis of Subsection 1.1 that whether a motion is bound or not 
depends both on the shape of the graph of the potential energy function and upon the 
value of the constant EF. Here the effective potential energy is given by Equation (11), 
and has the graph shown in Figure 7. The main features of the function U*® illustrated 
here are as follows. 


(i) When r is small, the term L?/(2mr?) dominates the term —GmM/r, so U®®(r) 
becomes large and positive as r approaches zero. 


(ii) When r is large, both of the terms in U*#(r) tend towards zero, so that the limit of 
U*"(r) for large values of r is zero. 


(iii) For a range of intermediate values of r, the term —GmM/r dominates the term 
L? /(2mr?), causing U°#(r) to be negative with a minimum value at some point 
rer". 


Exercise 6 


Find the value r = r* for which U*® has a minimum, and calculate the corresponding 
minimum value US, = U*#(r*). 

[Solution on page 50 | 

Figure 8 shows additions to the graph of U*® relating to two values of the total Bex 
mechanical energy E, which correspond to a bound and an unbound motion. A particle 
approaching the origin O with some positive energy, E¥>°%"4, reaches a closest point 
to O at r =a, following which it is always moving away towards larger values of r. On 
the other hand, a particle with some negative energy, E°°'"¢, undergoes a periodic 
oscillation of the radial coordinate between a minimum value r = b and a maximum 
value r = c. (In Subsection 1.1 these would have been bound and unbound motions 
along a line. Here they correspond to bound and unbound orbits.) 


Figure 8 


Exercise 7 
Describe the orbit corresponding to the total mechanical energy E = Lee 


[Solution on page 50 | 
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The following summary details what information has been obtained so far about orbital 
motion under the influence of gravity. 


Gravitational orbits 


Consider the planar gravitational orbit of a particle, with mass m and polar 
coordinates [r, 6], due to a spherically symmetric body of much larger mass 
M, situated at the origin. 


(i) The angular momentum L, of the particle is constant, where 


L, = mr76. (6) 


(This is a property of any central force.) 


(ii) The total mechanical energy E of the particle is constant. It is given by 


(a) — B= 4m(7? +1767) + U(r), (9) 


where U(r) = -GmM/r is the (actual) potential energy of the 


particle; 
(b) E = 3mr +U(r), 
where 
‘ ig GmM 
es Qmr2—r ) 


is the effective potential energy. 


(iii) The orbit is unbound if E > 0, and bound if Uct < E <0, where Ut 
is the minimum value of U*"(r). 


(iv) The equation of motion for r(t) is 


i GmM 
ea (8) 


r2 


mr = 


mr3 


The properties of gravitational orbits arising from the constancy of angular momentum 
and of total mechanical energy can be applied to obtain useful information in specific 
cases, as illustrated by the exercises below. 


Exercise 8 


A particle of mass m is launched with speed v at right angles to the line which joins it to a 
spherically symmetric body of much larger mass M. At the instant of launch, the distance 
between m and the centre of M is R. 


(i) Find the angular momentum of the particle for all subsequent times. 
(ii) Find the total mechanical energy of the particle for all subsequent times. 


(iii) Show that the motion will be bound provided that v* < 2MG/R. 


Exercise 9 


A particle of mass m is launched as described in Exercise 8. Its speed of launch v and distance 
of launch R are related by the condition v? = MG/(2R). 


(i) Find the total mechanical energy, and show that the resulting orbit is bound. 


(ii) Find ryin and rmax, the respective smallest and largest distances of the particle from the 
centre of the attracting mass M. (Hint: From the solution to Exercise 8, the total 


mechanical energy is E = >mv- — GmM/R, but this is also expressible as 

E = imr? + U*"(r) for all values of r and 7 on the particle’s trajectory. At the turning 
points, when r = rypin or r = Tmax, we have * = 0. Hence the equation for the turning 
points is E = U*#(r).) 


[Solutions on page 50 | 


Here’?,;;;, 800 Prone 
correspond to the values b 
and c of Figure 8. 
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3.3 The equation of gravitational orbits | 


We now revert to consideration of the radial equation of gravitational motion, 


£2 GmM 
(8) 


As mentioned earlier, we seek to obtain from this equation a relation of the form 

r =1r(@) to describe gravitational orbits. To do so, we apply the chain rule of 
differentiation to convert Equation (8), which specifies a function r(t), to an equation 
for r(@). For the first derivative, 7, we have 


en dr dear -dr 


"gt 6 a0 ae 


From Equation (6), 6 can be replaced by L,/(mr?), where L, is a constant. The 
previous equation can then be expressed as 


mr = 
mre r2 


ag 
~ mr? dé 
or, in a form which will shortly be convenient, as 
ee ee Se 


The equation of motion contains 7, which may be written in terms of a derivative with 


respect to @ by a further application of the chain rule. Thus, we obtain 
| i 2 oe. oF 
r= 0—r 


ae ep do 


which, after another use of the relation 6 = L,/(mr?) and substitution for * from 
Equation (12), becomes 


ee 
ee m?r2 dO? \ r ]- 


Equation (8) can now be written as 


se (2) - L? GmM 


mr? dé? \r} — mr? eee 
Defining 
1 
(9) = (8) (13) 
gives, after some rearrangement, 
d*u m?MG 
ee pe (14) 


Equations (13) and (14) contain the information necessary to deduce the orbit r = r(6). 
You have met equations like Equation (14) many times. It is a second-order, 
constant-coefficient, inhomogeneous differential equation, whose general solution has 
two arbitrary constants. One way of writing the solution is 
m?MG 
ee 
The two arbitrary constants in this solution are the amplitude B and the phase 6p. (We 
assume that B is non-negative. This represents no loss of generality, since the sign of 
the term Bcos(@ — 69) may be changed by adding 7 to the arbitrary phase 6.) The 
orbital equation is obtained from Equation (15) on using Equation (13) to replace u by 
1/r. After multiplying through by 
L? 
ene Se 16 
m?MG soe 
and putting e = BI, the orbital equation can be written as 


l 
oe 1 + ecos(@ — 69). (17) 


u(@) = Bcos(@ — 09) + (15) 


For example, in Unit 7 
Subsection 2.4 a similar 
equation was derived when 
modelling the motion of a 
particle suspended from a 
perfect spring in a uniform 
gravitational field. 


Here B is a non-negative but 
otherwise arbitrary constant, 
so the same is true of e. 
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Equations (16) and (17), together with the knowledge that L, (or, equivalently, 1) and 
e are constants, provide the required mathematical description of gravitational orbits 
for the model outlined at the beginning of this section. The remainder of the section 
will be spent in examining the implications of these results. 


Note first that the introduction of the symbols / and e has permitted us to write the 
orbital equation in a form which is identical to Equation (8) of Section 2. Provided that 
e <1, therefore, Equation (17) represents an ellipse with eccentricity e and semi-latus The terms ‘eccentricity’ and 


rectum I. ‘semi-latus rectum’ are used 
to describe e and l whether 
Exercise 10 or not the orbit is an ellipse. 


Equation (17) describes a general gravitational orbit. 


(i) By considering the possible values of the right-hand side of this equation, show that the 
distance of closest approach of the particle of mass m to the origin occurs at the angle 
? = 00, and has the value 


l 
l+e 
(ii) Show that, when e < 1, the maximum distance of the particle of mass m from the origin 
occurs at the angle 6 = 6) + 7, and has the value 
l 
eas 


(iii) Show that when e > 1 there is no maximum value of r. 


nin = 


Tmax = 


iv) From parts (i)—(iii , give conditions on the eccentricity e which characterize bound and 
gS y 
unbound orbits. 


[Solution on page 51 | 


From the results of Exercise 10 we can sketch the possible types of orbit that 
Equation (17) represents. This is done in Figure 9 below. When e < 1, the orbits are 
closed (ellipses) with closest and furthest distances rmin and Tmax respectively. When 
e > 1, the orbits are unbound (hyperbolas or parabolas), and only the distance of 


closest approach rmin is defined. In either case, the orbit is symmetrical about the line 
0 — Ao. 


Figure 9 


It is convenient to consider the gravitational orbits with respect to a set of axes which 
are rotated anticlockwise by an angle 0) from those of Figure 9. This is equivalent to 
choosing 99 = 0 in the orbital equation, or to choosing the x-axis to coincide with an 
orbital symmetry, and represents no loss of generality. In place of Equation (17) we 
then have 


l 
~ =1+ecos8, (18) 


where the semi-latus rectum / is given as before by Equation (16). The corresponding 
pictures for the orbits are shown in Figure 10. 
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min 


(aje<l (Byes 


Figure 10 


Equation of gravitational orbits 


The Newtonian model for the gravitational orbit of a body of mass m around 
a spherically symmetric body of much larger mass M located at the origin 
leads to the orbital equation | 


l 
-=1+ecos8, (18) 
T 


where the eccentricity e is non-negative and the semi-latus rectum | is given 
by 
| - 
e Se 16 
m?MG /» 
Here L, is the constant angular momentum magnitude of the orbiting body. 
If 0 < e < 1 then Equation (18) describes an ellipse. 


Values for the constants / and e in Equation (18) can be calculated in a specific case 
from experimental observations, and this can be done in various ways. For example, 
suppose that the distance of closest approach rin and the corresponding angular 
velocity 6 have been measured, and that the value of the central mass M is known. 
Then a value for L,/m can be obtained from Equation (6). This leads in turn to a 
value for / from Equation (16) and a value for e from the result of Exercise 10(i). 


The energy analysis carried out in Subsection 3.2 predicted that gravitational motion This analysis was based on 
will be bound if the total mechanical energy E satisfies the inequality Uc! < E < 0, Equations (10) and (11), 


and unbound if E > 0. It is to be expected that these predictions are related to those together with Figure 8. 


obtained in terms of the eccentricity e in Exercise 10, namely, that an orbit is bound if 
e <1 and unbound if e > 1. We consider this relationship in the example below. 


Example 1 


Show that the total mechanical energy can be written as 


mMG 
ES *— 1). 
Hence verify that the condition EF > 0 for unbound orbits is equivalent to e > 1, and 


that the condition E < 0 for bound orbits is equivalent to e < 1 (given that e > 0). 


Solution 
From Equation (10), E can be written as 
L? GmM 


E=1mr? + 
2 2mr2 r 
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Using Equation (12) for 7, this becomes 
Pe ee ee 
Baim (i = 
m? |d0\r 2mr? r 


From Equation (16) for 1, this is 


e=He (Te ()] +(4)'-2(} 


But l/r is given in terms of 6 by Equation (18), so that 


M 
= er {e* sin? 6 + (1 + ecos)? — 2(1 + ecos6)}. 
On simplifying, this gives | , 
mMG 


21 


So the condition E > 0 for unbound orbits corresponds to e > 1, and the condition 
E <0 for bound orbits corresponds to e < 1. (Although this outcome is expected, it 
represents a useful cross-check on the results achieved.) O 


Exercise 11 
Describe the orbit given by Equation (18) when the eccentricity e is zero. What are rmax and 
Tmin 1n this case, and what is the total mechanical energy? 


Exercise 12 


A satellite of mass m follows a circular orbit with constant angular velocity and radius d 
around a spherically symmetric planet having a much larger mass M. 


(i) Show, from first principles (using Equation (13) of Section 1), that the satellite’s 
acceleration lies along the radius of its orbit and has magnitude 


|é| = dé. 
(ii) Show too that its orbital angular momentum has magnitude 
C= md?6. 


(iii) By equating m# to the gravitational attraction on the satellite (Newton’s second law) 
show that the value of the orbital radius agrees with that found in Exercise 11. 


(iv) By adding the kinetic energy of the satellite to its gravitational potential energy, verify 
that the value of the total mechanical energy agrees with that found in Exercise 11. 


[Solutions on page 51 | 


3.4 Planets and comets 


In this brief subsection we shall verify the statement made at the beginning of this 
section to the effect that Kepler’s laws of planetary motion are consequences of our 
model based on Newton’s laws. Following this there is a short discussion about the 
orbits of comets. 


Kepler’s first law says that each planetary orbit is an ellipse, with the Sun at one focus. 
You showed in Exercise 3 that a planet’s orbit lies in a plane because gravity isa 
central force. According to the model, a gravitational orbit within this plane has an 
equation of the form 


] | 
= 1+ ecos8, 3 (18) 


where the semi-latus rectum / and the eccentricity e are positive constants for a 
particular planet. 


Planetary orbits are clearly bound, and you established in Exercise 10 that e < lisa 


necessary and sufficient condition for bound motion. But for e < 1, Equation (18) is the 


equation developed in Subsection 2.1 to describe an ellipse with one focus at the origin. 
In our Newtonian model of planetary motion, the Sun is located at this focus. Thus 
Kepler’s first law is a consequence of the model. 


Equation (7) of Section 2 
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Kepler’s second law says that the line joining a planet to the Sun sweeps out equal areas 
in equal times. We showed in Subsection 3.1 that this is equivalent to the equation 


5 (076) =0. (3) 


which expresses the constancy of the planet’s angular momentum. This is again a 
consequence of gravity being a central force. 


Kepler’s third law is a relation between the period of a (bound) orbit and its lateral 
_ extent. It says that if 7’ is the period and a is the semi-major axis of the elliptical orbit, 
then | 

To = he (19) 


where k is a constant not depending on T or on a. In Exercise 13 below you are asked 
to show from basic principles that Equation (19) holds for the special case of a circular 
orbit, where the eccentricity e is zero. In this case the semi-major and semi-minor axes — 
both equal the radius of the circle. In Exercise 14 you are then asked to establish 
Kepler’s third law for any elliptical orbit. 


Exercise 13 


(i) Show that the magnitude of the gravitational force between the Moon (of mass m) and This exercise refers to 
the Earth (of mass M) is Newton’s quotation given at 
GmM : the beginning of 
oe Subsection 2.2. 


where d is the separation between their centres. 


(ii) Assume that the Moon is moving uniformly in a circle around the Earth, so that (from 
Exercise 12(i)) the magnitude of its acceleration is |r| = d0?, where 6 is its constant 
angular velocity. Show that 

4nd 
Ir | = T2 ? 
where T' is the period of the Moon’s orbit. 


(iii) Show that if air resistance is ignored then the acceleration of an apple dropped close to 
the surface of the Earth is 
aS 
= 
where R is the Earth’s radius. 


(iv) From parts (i)—(iii) show that, according to the Newtonian model of gravitation, 


2 2 An? 
T= he where = gh 
(v) It was known by Newton that d ~ 60R, g ~ 9.8ms~? and T ~ 28 days, but it is now 
believed that he used for R, the Earth’s radius, a value of 3500 ‘Italian miles’, where one 
Italian mile is 5000 feet (rather than 5280 feet). This would make R ~ 5.3 x 10° m. Using 
these values (converted to appropriate SI units where necessary), do you find that the 
left- and right-hand sides of the predicted relation in part (iv) agree ‘pretty nearly’? (A 


more accurate value for FR is 6.4 x 10° m.) 


Exercise 14 


Consider an elliptical orbit with semi-major axis a and semi-minor axis b. Then it is known 
from Exercise 2 of Section 2 that the area of the ellipse is tab. Also, the semi-latus rectum is 
b? ee 
ae a m2MG 
(from Equation (5) of Section 2 and Equation (16) of this section), and the rate at which area 
is swept out is | | 


ee 
A=65r — 


Ls 
2m 


(from Equation (9) of Section 2 and Equation (6) of this section). Here L. is the (constant) 
orbital angular momentum. 


(i) From these facts, show that the orbital period, T, and the semi-major axis, a, are related 
by Kepler’s third law, in the form 

An* 

7 MG 

(ii) Show that, for a circular orbit around the Earth, this agrees with the result of 
Exercise 13(iv). 


T? = ka’, where k= 


[Solutions on page 51 | 
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Kepler’s laws were formulated to describe the bound orbits of planets, so there is no 
reason to suppose that they should have significance for unbound orbits. The first and 
third laws can by their very nature apply only to the bound case. However, Kepler’s 
second law applies even for unbound orbits since, as you have seen, it depends only 
upon the fact that gravity is a central force. 


The Newtonian model which led to Equation (18) for gravitational orbits makes no 
assumption as to whether or not the motion is bound (with eccentricity e < 1) or 
unbound (with e > 1). As shown in Example 1, these conditions for the eccentricity are 
related to similar conditions on the total mechanical energy E of an orbit. If EF < 0 for 
a body moving under the gravitational influence of the Sun, then the body has 
insufficient energy to ‘escape’ from the Sun’s gravitational field and is confined to 
bound and periodic motion. Bodies with energy E > 0 may approach the Sun once, but 
are not ‘captured’ by it. These unbound orbits are called parabolic when e = 1 (E = 0) 
and hyperbolic when e > 1(E > 0). 


At any particular time there are hundreds of bodies in the solar system moving in nearly 
or actually unbound orbits under the Sun’s gravity. These are the comets. Precise 
details of cometary orbits are not easily obtained, but it appears that most comets are 
in weakly bound orbits, with values of eccentricity e just less than unity. These comets 


therefore have highly eccentric, or flattened, orbits, as illustrated in Figure 11. You investigated the 
connection between the 
eccentricity and shape of an 
ellipse in Exercise 4 of 
Section 2. 


Figure 11 


As a first approximation, the orbits of these bound comets obey Kepler’s laws, but as 

they follow their elongated trajectories through the solar system they occasionally pass 

close to a planet. In these circumstances the gravitational effect of this third body can 

alter their orbits, even to the extent of raising their energies to positive values. Owing 

to its relatively large mass, Jupiter is a particular offender in this regard. The most 

famous bound comet is Halley’s comet, which has a period of about 76 years and last 

passed close to the Sun in early 1986. Its orbit has eccentricity e ~ 0.97 anda 

semi-major axis of about 18 times the distance from the Earth to the Sun. Comparable 

data for the planets, the Moon and Pluto are given in the table below. For each body except the 
Earth, the values given in the 
table for mass, period, 
semi-major axis and mean 


diameter are multiples of the 
corresponding values for the 


San 332 948.0 ie oe 
Moon 0.012 771 
Mercury 0.055 0.241 0.387 0.391 
Venus 0.815 0.615 0.723 0.951 


Period of Semi-major Eccentricity Mean 
revolution axis of orbit of orbit diameter 
around Sun 


Earth 5.974 x 1074kg | 365.256 days | 149.6 x 10®°km 12 735km 
Mars 0.107 1.88 1.524 0.531 
Jupiter 317.9 11.86 5.203 10.98 
Saturn 95.1 29.46 9.539 9.09 
Uranus 14.6 84.00 19.182 4.00 
Neptune 17.2 164.79 30.058 3.78 
Pluto 0:17 247.70 39.44 0.22 


At any time there are comets in the solar system that are pursuing unbound 
trajectories. These bodies come in from distant places, pass partly around the Sun to a 
distance of closest approach and then recede, never to return to our solar system. 
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Summary of Section 3 


1. 


A central force F, acting from the origin on a particle with position vector r, has 
the form 


F(r) = g(r)r, 


where g(r) = g(x,y, z) is a scalar function of position. If the particle has mass m 
and is not acted upon by any other force, then it has constant angular momentum 
L=r x mr. If the unit vector k is chosen in the direction of L, then 


(i) the motion of the particle is confined to the (z, y)-plane; 
(ii) L=L zk, where L, = mr?60 is constant. 


The gravitational force due to a spherically symmetric body whose centre is fixed 
at the origin is a central force. In this case, the constancy of L, is equivalent to 
Kepler’s second law. 


An isotropic central force has the form 
F(r) = f(r)e,, 
where r = ,\/z? + y? + z? is the distance from the origin, and e, = r/r. Fora 


particle of mass m acted upon by such a force, the radial equation of motion is 
ep 


mr3 


mr = f(r) + 


which may be integrated (analytically or numerically) to give r(t). The 
gravitational force is isotropic, with f(r) = —-GmM/r?, where M is the central 
mass and G is the gravitational constant. (It is assumed that M is much larger 
than m, in order that the central mass should remain fixed at the origin.) 


For the particle moving under gravity as described above, the total mechanical 
energy F is constant. It may be expressed as 


(i) E = 1m(r? +776?) + U(r), 

where U(r) = —GmM_/r is the (actual) potential energy of the particle; 
(ii) E=imr?+U(r), 

where 


L2 GmM 


2mr? r 


u(r) = 


is the effective potential energy. 


Orbits are unbound if & > 0, and bound if Vel < E <0, where US#. is the 
minimum value of U*"(r), 


The equation of gravitational orbit for the particle described above is 


l 
—=1+ecos8, 
2 
where the semi-latus rectum / is given by 
L? 
i= ca 
~ m*MG 


If 0 < e < 1 then the orbit is an ellipse with eccentricity e, showing that Kepler’s 
first law is satisfied. If e > 1 then the motion is unbound. 


In terms of the eccentricity and the semi-latus rectum, the total mechanical energy 
is 
GmM 


= 
21 


(e=1). 


Kepler’s third law is also a consequence of this Newtonian model, since the period 
T of an elliptical orbit and its semi-major axis a are related by the equation 


Arr? 


T? =ka*®, where ade v7 
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4 Central forces and conservation laws 
(Television Section) 


4.1 Isotropic central forces 


The television programme considers some general statements that can be made about 
the motion of a particle in three dimensions which is acted upon by a central force and, 
more especially, by an isotropic central force. When the force is central, then the 
angular momentum of the particle is constant and its motion lies in a plane. When in 
addition the force is isotropic, one can say further that there is a potential energy 
function for the particle, so that it has a constant total mechanical energy. 


The first of these results was established in Subsection 3.1. The second was applied for 
the particular case of gravity in Subsection 3.2, but holds more generally. This more 
general version is a consequence of Result (i) below. 


(i) Ifa scalar function ¢(r) depends upon position only through the radial 
distance |r| = r = (x* + y? + z?)!/?, so that o(r) = d(r), then the 
gradient of ¢ can be written as 


grad ¢@ = a 
dr 


where e, = r/r is the unit vector in the outward radial direction. 


(ii) Hence if F is an isotropic central force, that is 
Fir) = f(r)e, 


for some function f, then F is conservative, with a potential energy 
function U(r) such that 


dU 
DS ee 


The motion of a particle in this force field is planar, and has constant 
total mechanical energy 


E = 3m(r? + r267) + U(r). 


Result (ii) follows from Result (i) on referring to the definition of a conservative force 
(Equation (3) of Section 1) and to the formula for E derived in Exercise 13 of 
Section 1. (The motion is planar because the force is central, as in Subsection #1.) 


The proof of Result (i) is achieved by applying the definition of the gradient and the 
chain rule. Thus we have 


Since ¢(r) = @(r), this can be written as 


dp Or, ddoOr, ddoOr gO-1 Or. dr. Or 
dg= — i+ ——j+ ——k= — {| —i+ —j+ — 
eee dr Ox dr ay? dr a dr aye - 


But from r = (x? + y? + 27)!/?, we obtain 


Le ee te eee ee oye x 
Pe + y + 2°) = 31% Sas + 2°) ee 


Similarly, we find that 
OF. ¥ oF 3 


es d 2 
Oy or ee a 
so the formula for the gradient becomes 
d 
grad ¢ = a (=i + = — =k) 
ar Ne a r 


Lt ee dpr do 
a k = — - = - 
sae aes ) dr r Pe 
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This completes the proof of Result (i). 


Now watch the television programme. 


4.2 Programme synopsis 


Tom Smith introduces central forces. As an example, he demonstrates an almost 
frictionless puck, connected by a spring to a fixed origin but otherwise free to move on 
a horizontal surface. A central force acting on a particle is a force which is directed 
either parallel or antiparallel to its position vector r. Such a force can be written 
mathematically as 


F(r) = g(r)r, 


where g(r) = g(x,y, z) is a function of position. It is shown that the angular 
momentum of the particle, L = r x mr, is conserved (constant in time) if the particle is 
acted upon only by a central force. In this case the particle moves in a plane, which 
may be taken as the (2, y)-plane. With this choice of axes, the constunt vector L is 
directed along the z-axis. The fact that a force is central thus allows something to be 
said about the corresponding particle motion, even when it is not possible to solve its 
equations of motion analytically. 


Mike Crampin takes up the story by pointing out that, since L is constant (L = 0), 
both its direction and magnitude are constant. The motion is in the (x, y)-plane and, 
using plane polar coordinates, the angular momentum is 


| Pea ee 


where L, = mr. 


Since the magnitude L, is constant, this relation can be used to eliminate either r or 6 
from the equations of motion. The constancy of r2@ is equivalent to Kepler’s second 
law, which says (when the force involved is that of gravity) that equal areas are swept 
out in equal times. So this law of Kepler, together with that part of his first law which 
implies that planetary orbits are planar, follows from the fact that gravity is a central 
force. 


At the next stage, zsotropic central forces are introduced. These are central forces 
whose magnitude depends only upon the distance |r| from the origin. Such forces have 
the form 


F(r) = f(rjer, 
where r = |r| and e, is the unit outward radial vector. This is a special case of the 


general central force, g(r)r. 


At this point the second result proved in Subsection 4.1 is applied. ‘This says that if a 
force is isotropic as well as central then the total mechanical energy of the particle is 
conserved as well as its angular momentum, for then there exists a potential energy 
function U(r) such that 


d 
F = = gradU = Ps 
dr 


Using plane polar coordinates to describe the motion, the total mechanical energy is 
E = im(?? +1767) + U(r). 


But i. = mr26 is constant too, so # can be written as 


This equation expresses the energy solely in terms of r and 7, resembling the energy 
equation for a one-dimensional motion, so the radial motion can be considered 
separately. To make use of the one-dimensional analogy, it is convenient to define an 
effective potential energy function, U*", by 
L? 
pet (r) oe z 


2mr2 


+ U(r). 
Then the constant total mechanical energy is 


= imi + U(r). 


TV30 


Equation (1) of Section 3 


See Subsection 3.1 for details 
of this analysis. It is assumed 
here that L 40. You 
investigated the case L = 0 
in Exercise 3(iii) of Section 3. 


See Exercise 4 of Section 3. 


Equation (4) of Section 3 


In fact, the potential energy 
function is 


U(r) = - ; f(s) ds, 


where ro is any non-negative 
constant. 


The constancy of E (FE = 0) 
was shown in Subsection 1.2, 
and its expression in plane 
polar coordinates was derived 
in Exercise 13 of Section 1. 
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Looking at the graph of the function Ue", and using the constancy of EF, allows us to 
see whether there are any turning points for the radial motion. At these points 7 
vanishes, and the radial motion reverses direction. The details depend upon the precise 
form of the potential energy function U and on the value of E. If there are two turning 
points then these are the closest and furthest distances from the origin reached by the 
particle at the given value of energy E. If there is one turning point only then the 
motion is unbound, and the turning point is the distance of closest approach. 


If the potential energy looks like that in Figure 1(a), then U®# is given by the graph in The potential energy given in 
Figure 1(b) and there are always two turning points, for any energy /. But for a Figure 1(a) is that for the 
potential energy like that of gravity, for which the graph of U is that in Figure 2(a), 
there will be bound or unbound motion, depending on the energy level. This is 
illustrated in Figure 2(b). 


puck attached to a perfect 
spring, with stiffness k and 
natural length ro. 


Figure 1 


Lett 


Eenbound 


For gravity, the potential 
energy is 
eas GmM | 
r 
as in Equation (19) of 
Section 2. 


(a) (b) 
Figure 2 


The puck connected to a spring is reconsidered in the light of this analysis. Since 
friction is negligible the force is central, and the effective potential energy looks like 
that of Figure 1(b). From this graph, inner and outer bounding radii (radial turning 
points) are predicted for the motion. Viewing the puck’s motion from an overhead 
camera verifies this prediction. 


In summary, isotropic central forces conserve both angular momentum and total 
mechanical energy. Without solving the equations of motion, the turning values for r 
can be found from knowledge of the potential energy function U (r) and total 
mechanical energy EL. 
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Two further physical examples are considered. The first is gravity, which is the main 
topic of this unit. The gravitational potential energy U and its corresponding effective 
potential energy U°" are the functions shown in Figures 2(a) and 2(b) respectively. So 
in this case there may be two turning points, as for elliptical orbits, or only one turning 
point, as for unbound hyperbolic orbits. 


The second example is the mutual repulsion between two electrostatic charges of equal 
sign, either both positive or both negative. Here the potential energy U has the same 
general shape as the gravitational case (Figure 2(a)), but with opposite sign. The 
functions U and U°" for the electrostatic force are shown in Figures 3(a) and 3(b) 
Lelow. It is clear that there are no bound states in this case. The potential U for both 
examples can be written in the form U = y/r, where pu < 0 for gravity and pu > 0 for 
electrostatic repulsion. 


U Uet 


Founbound 


Figure!3 


The physicist Ernest Rutherford carried out an experiment in the first decade of the 
twentieth century to demonstrate that atoms consist of relatively massive, extremely 
small nuclei bearing a positive electrostatic charge, surrounded by light electrons 
carrying a compensating negative charge. The experiment consisted of firing the lighter 
nuclei of helium atoms, with some energy EF, at other heavier atoms and measuring the 
(perforce) unbound ‘scattering’ orbits of the helium nuclei. Rutherford achieved good 
agreement with his experimental data by applying the classical mechanical theory of 
orbits. Surprisingly, the modern theory of quantum mechanics gives similar theoretical 
predictions, for subtle reasons. 


Summary of Section 4 


1. Ifascalar function ¢(r) depends upon position only through the radial distance 
lr] =r = (27+ y? + 27)!/?, so that d(r) = d(r), then the gradient of ¢ can be 
written as 

d 
grad ¢@ = Lae 
dr 
where e, = r/r is the unit vector in the outward radial direction. 
2. Hence if F is an isotropic central force, that is, 
Fir) = f(rje, 


for some function f, then F is conservative, with a potential energy function U(r) 
such that 


The motion of a particle in this force field is planar, and has constant total 
mechanical energy 


E = im(7? +1776?) + U(r). 
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5 End of unit exercises 


Exercise 1 


For a gravitational orbit, it was shown in Example 1 of Section 3 that the total mechanical 
energy is 


2 
e'-1 

2] ( ), 

where m is the mass of the orbiting particle, M is the mass of the central body, G is the 

gravitational constant, and e is the eccentricity of the orbit. 

(i) For a bound orbit, show that 


mMG 
7 
where a is the semi-major axis of the elliptical orbit. 


B= 


(ii) For any orbit, show that 


where rmin is the distance of closest approach of the orbiting particle to the origin. 


Exercise 2 


It is proposed to put a geostationary satellite into a circular orbit above the Earth’s equator, 


such that the period of the satellite’s orbit matches exactly the period of the Earth’s rotation. 


Show that the distance d of the satellite from the Earth’s centre is given by 
2p = gR°T* 
An2 ° 
where R is the Earth’s radius and T is the period of its rotation. Taking g = 9.81ms ? and 


R=6.4 x 10° m, calculate the corresponding height of the satellite above the surface of the 
Earth. 


Exercise 3 


A particle of mass m is launched with speed v at right angles to the line which joins it to a 
spherically symmetric body of much larger mass M. At the instant of launch, the distance 
between m and the centre of M is R. You showed in Exercise 8 of Section 3 that, for the 
subsequent motion, the angular momentum has magnitude 


Li, = nek 
the total mechanical energy is 
de imo" = ae 
and the motion is bound if 
s 2MG 
ie 


(i) For bound motion, show that the semi-major axis of the orbit is 


2 v" = 
=e (G 2 aa) 


(ii) For bound motion, find an expression for the eccentricity e in terms of v, R, M and G. 
(iii) Show that the condition for circular motion is 
25: MG 
= 


[Solutions on page 52 | 


This situation was considered 
previously in Unit 28 
Subsection 4.2. 


MST204 30 Solutions 


AT 


Appendix: Solutions to the exercises 


Solutions to the exercises in Section 1 


1. (i) Integration of dU/dx = —k gives U(x) = —kx +C. 
Here C’ is an arbitrary constant which may be chosen equal 
to zero for simplicity, leaving U(x) = —kz. 
(ii) Integration of dU/dz = kx gives U(x) = $kx* +C, 
where C’ is an arbitrary constant. Choosing C’ = 0 produces 
Uta) = tke. 
(iii) When U(x) = k/z, the force is 

F(x) = -dU/dx =k/x* (xz > 0). 


2. Differentiation of 
a(t) = 29 coswt + -° sinwt 
Ww 
gives 
z(t) = —waxo sinwt + vo cos wt. 
(Note here that +(0) = zo and «(0) = vo, as claimed.) 


Squaring each of these equations gives 
2 


2 i. 2 2X0 Vo ' vo . 2 
tf =, COs wl cos wt sin wt + —z Sin’ wt 
Ww 
and 
“2 ea ace . 2 2 
Lr =w x£osin wt — 2wxrovo cos wt sin wt + vg cos” wt. 


The total mechanical energy is 
Et = ima” + Lkax* +C 
= ima” + imw’2? +C. 
On substituting for x? and <? from above and simplifying, 
this becomes 


E(t) = $mw 29 + 3 liu, + C = BO): 


Hence, by direct substitution, it has been shown that F is 
constant. 


At various times the particle reaches its turning points, 
located at x7 say. Since the speed vanishes at these points, 
the energy there can be expressed as 


Ee= ikap+C = 


Hence, since F is — oe the motion, we have 


mw t+C= tm” x + imvg +C, 


from which the locations of the turning points are deduced 
to be 


Poe ae & 


mw “75 PS. 


re + ve /w?. 


3. Since the particle is released from rest, its initial velocity 
x(0) is zero, so that its total mechanical energy is then 


E = imz*(0) + U(a) = U(a). 


The energy is constant throughout. At 7 = (@ it can be 
written as 


E = 5mv (6) + U(8), 
where v(x) = «. Equating the two expressions for EF’ gives 


pmv (3) + U(B) = U(a), 
so the speed at x = (@ is 


u(B) = 4/ — (Ula) — U(A)). 
Similarly, the speed at x = 7¥ is 
= (U(a) - U(>)). 


For x > 6 the potential energy function U(x) is zero, so the 
speed is then constant, with value 


From Figure 6 of Section 1, or by considering these formulas 
for the speed and remembering that U() < 0,U(y) > 0, we 
have 


v(B) > v(6) > v(7). 


4. The particle is projected from a point to the right of 
x = 6, where the potential energy vanishes (see Figure 6 of 
Section 1). Since its initial velocity is 7(0), the energy is 


E = imz’(0), 


from which we have 

SSS bala 

m 

The negative sign here arises because the particle is 
projected to the left. After projection it slows down toa 
minimum speed at x = y and then accelerates to an overall 
maximum speed at x = GZ. After reaching this point it slows 
to a halt at r = a. Thereafter it retraces its path, with a 
maximum speed at x = 7 and minimum speed at x = ¥. 
After reaching x = 6 it moves out to infinity with velocity 
equal to its initial speed, that is, 


2B 

=e 

5. Using the definition of the arenes the vector equation 
F = — grad U, 


with F = —mgj, is equivalent to Lie three scalar equations 
oe He 
= ee Sree Sa 


By inspection, these equations have the solution 
U=mgy+C, 


where C is any constant. Thus the uniform gravitational field 
is conservative, with the potential energy function just found. 


6. We seek 
oa ne 
F = — gradU = ———i —- —j -—- —k 
se Ox By? ee 
where 
U =k/r=k(a? +9" +27) 
Now we have 


i792 


a = —tk(x* pa + Pies ce x 22 = —kar* 
x 
Similarly, it can be shown that 
a =-—kyr* and — = —kzr*, 
so that 
k; k 
= 3 (ai +yj+ zk) = = 


7. It is given that 
mr=-—gradU+G and G-r=0. (1) 


Differentiating E = mr + U and applying the chain rule 


leads, as in the text, to the equation 


E=r-(mt+gradU). 
Using each of Equations (1) in turn, this becomes 
BE =4G =. 
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8. As explained in the text prior to the exercise, the total 
mechanical energy FE is constant, where 


a imu” +mgytC. 
Initially, the particle’s speed is zero, since it is released from 
rest, and its height is y = h, so that E = mgh+ C. When it 
reaches the origin we have a> = foes u = uo, so the energy 
can also be written as EF = sMup + C. Hence 

Imus +C = mgh+C, : 
giving the particle’s speed at the origin as 


to = +7 2gh. 


(Note that this outcome is independent of the slope of the 
plane.) 


9. Since i,j form an orthogonal pair of unit vectors, we have 
i<j =D. ane +4 = 5-3-4, 
Then Equations (9) of Section 1 give 
e, + e€9 = (cos#i + sin j) - (— sin8i + cos@j) 
= — cos @sin 0i-+i+ (cos” 6 — sin? 0)i-j 
+ sin@cos6j-j 
so); 
Similarly, we obtain 


e, +e, = (cos i+ sin 8 j) - (cos6i+ sin 6 j) 
= cos’ +sin? 9 =1 
and 
eg - eg = (—sin#i+ cos@j) - (—sin8i + cos @j) 


— sin? 0 + cos? 6 = 1. 


10. Since the particle is constrained to move at a fixed radius 
r = R, we have r = 7 = 0. Then, in plane polar coordinates, 
—R6’e, + Roe. 

(These are the results which you obtained in Exercise 7 of 
Unit 28 Section 4, and which were derived again in 


Subsection 4.2 of that unit.) Using Equations (9) of 
Section 1 to substitute for e,,eg in terms of i J gives 


r= —ROsin0i+ RO cos 0j, 
r=-R (6° cos 6 + 6sin6) i+ R(@cosé — 6° sin 8) j J. 


r= R6ep and *+= 


11. The polar angle 9 is fixed here at the value 69. Hence 
6=6=0. Then, in plane polar coordinates, 

r= 76, 2h = 7e,. 
In terms of the Cartesian unit vectors i and j, these are 

r = 7rcos69i+7sin 4 j, 

r = Fcos 091+ 7'sin Oj. 


12. We have : 
r= at’ + bt, f= 2ot4 bh, # = 2a. 
and 8 ct”, @ = 2ct, O= Ze, 


so that the velocity and acceleration are given by 
r = (2at + b)e, + (at® + bt)(2ct)e, 
and ¥ = [2a — (at? + bt)(2ct)"Je, 
+ [2(2at + b)(2ct) + (at? + bt)(2c)Jes 


13. The velocity of the particle is 
r=re, + rbeo, 
and since e,,eg are orthogonal unit vectors (with e, -eg = 0 
and e, -e, = e9-eg = 1) the total mechanical energy is 
B= 4b + 


2 


= 5m (** +1776?) +U. 


14. (i) The force acting on the particle is F = —Te,. Since 
r is a linear function of t, we have 7 = 0. Hence Newton’s 
second law, in the form of Equations (14) of Section 1, gives 


—mré? =—T and ae (r°0) =f, 


t 
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These equations are equivalent to 
mré” = T, (2) 


d (2 
=, (r°8) =0. (3) 


(ii) From Equation (3), we know that r76 is constant. 
Hence we have 


r”(t)0(t) = r?(0)0(0). 


This is used to substitute for 6(t) in Equation (2), giving 


279\f 2 4/\f2 4/1\ p2 
Q , 
P= tae) r~(0)0(0) _ mr (0)0“(0) _ mr (0)A-(0) 
r?(t) r3(t) (r(0) — ut)? 
as required. 
(iii) The string will break when T = 500N. From the answer 
to part (ii), this occurs when 
4(n\ 92 
mr” (0)8* (0) ~ 500, 
a 
or when 


r = (mr*(0)6?(0)/500) /°. 


With m = 0.1kg, r(0) = 1m and 0(0) = 1rads~’, this gives 
r = (0.1/500)'/? ~ 0.058m. The time at which the string 
breaks is then 

t = (r(0) — r)/p = 60(1 — r) ~ 56.5s. 
(Note, from the form of the function for 7’, that for any finite 


breaking strength the string will break before the particle 
reaches the origin.) 


15. (i) The force is —k(r —1,)e, so that, from 
Equations (14) of Section 1, the equations of motion are 


m (7 - 70”) = —k(r — lo), 
2 
ae (r°0) = {j 


(ii) From the last equation, mr?6 must be constant. (This 
quantity is the angular momentum of the particle. 
Conservation of angular momentum occurred also in the 
previous exercise, and will be of great importance in the 
derivation of planetary orbits.) 


(iii) If U = $k(r — 1o)*, where r = (x? + y”)'/?, then we have 
OU Or 
oa per sk x A(r — lo) a 
= k(r — Io) x L(a? +y jes 
= k(r — Ip)ar™* 
Similarly, 
OU 
ay = k(r = lo)yr = 
so that (using Equation (8) of Section 1) 
OU. OU , 
peel <2 Se oe 
grad U = ae Dy 
= —k(r — lo)r~* (xi + yj) 
= —k(r —lIo)r-tr 
= —k(r aie lo )e,. 


Since this expression for — grad U is equal to the force F 
given in the statement of the exercise, we have shown that U 
is a potential energy function for F (and hence established 
that F is a conservative force). 


(iv) Suppose that the speed of the particle is denoted by v. 
From part (iii), the total mechanical energy 


E = imv* +U= mv * + Ek(r — Io)? 


is constant. When r = 2lo, the speed is given as v = wu, so 
that 


i= imu” LklG. 
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When r = lo the potential energy is zero, so that the speed v 


is then given by F = mv. Hence we have 
fo ee 
smu = smu’ + 5klo 
ce 
or v= (/u2+—2 
m 


for the speed of the puck when the spring has its natural 
length. 


Solutions to the exercises in Section 2 


1. When a = b, the equation of the ellipse can be written as 
XoLY =a 
This is the equation of a circle with radius a. 


2. (i) The area A is twice the area of the upper half of 
the ellipse, namely 


Anal Y(X)dX. 


a 


But from the definition of the ellipse (Equation (1) of 
Section 2), the upper half of the ellipse (where Y > 0) has 


1 — (X/a)?. Thus the area is 


A= 2» f i (var ae 


equation Y = b 


(ii) Let X =acosu. When X = —a, u=7, when X =a, 
u = 0, and dX = —asinudu, so that the integral becomes 


‘ 0 Tw 
A= -2ab | sin? udu = 2ab | sin? udu 


0 


= ab | (1 — cos 2u) du = zab. 
0 


3. <A circle corresponds to a = b, for which gq = 0. Hence, 
for a circle, the foci coincide at the centre. 


4. (i) When }=<a the eccentricity is e = 0. 
(ii) When b = $a, we have e = $/3 ~ 0.866. 


2 


(iii) When b = 4a, we have e = +V8 ~ 0.943. 


The larger the semi-major axis a is compared to the 
semi-minor axis b, the more ‘squashed’ is the shape of the 
ellipse. This corresponds to larger values of the eccentricity e. 


5. From Equation (9) of Section 2, we have 


A = 4R70. 
But also V = RO, so that 
A= iRV. 


Since both R and V are constants, A is constant. 


6. Refer to Figure 4 of Section 2. During the interval 

ty <t < te, the planet is closer to the origin than it is during 
the interval ts < t < t4. So to sweep through the same area 
(as it must do by Kepler’s second law), it moves faster 
during the first interval. 


7. If T isthe planet’s period and a is its semi-major axis, 
then Kepler’s third law states that 


T? = ka®, 
where k is a constant. If 7\,am are respectively the period 


and semi-major axis for Mars, and 7, ag are the 
corresponding parameters for the Earth, then we have 


2 3 3/2 
(=) = (=) oo (=) 
Tr aR an 


For the values given, 
Tm = 365.256 x (1.524)?/? ~ 687.187 days. 
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8. In terms of units, Equation (10) of Section 2 for the 
force of gravity gives 
2 
ps G x (mass) | 
(length)? 
where ‘=’ here means ‘has the same units as’. It follows that 
(length)* x force 
G = 
(mass)? 


But from Newton’s second law, F = ma, we have 
mass x length 


fe = 
orce Gime: 2 
so that 
ae (length)? _—mass x length (length)? 


(time)? =—s mass x (time)?’ 
Thus, in SI units, G has the units m°’ kg~'s~?. 


~ (mass)? 


9. The given vector is of the form a/|a|, where a = r; — rj. 
Any vector of this form is a unit vector, since its magnitude 
is |a|/|a| = 1. 


10. Equation (10) of Section 2 gives 


Gmim; 
Fij = ees GF = 7p 


as the force on m; due to m;. Hence, swapping the indices 7 
and j, the force on m; due to m; is 


Gm;m; Gmim; 
jt Ir; <= r;|3 (r r;) Ir; =e r;|° (r; ri), 
where we have used the fact that |r; —r;| = |r; — r;|. Thus 


F;; = —F,;:, and Newton’s third law is satisfied by the 
gravitational force. 


11. Equation (15) of Section 2 gives the gravitational 
potential energy function resulting from a system of particles 


mi,™M2,...,mn, acting upon another particle of mass m 
located at r. Here there is only a single particle in the 
system (N = 1), with mass m; = M and position vector 
r; = O. Hence the potential energy function is 
wae _GmM = _GmM 
Ir| r 


The corresponding gravitational force is 
F(r) = — grad U(r). : 
From Equation (13) of Section 2, with r; = r; = 0, this gives 
i = Geller 
ee ee’ 
which could also have been obtained directly from 
Equation (12) of Section 2. Putting e, = r/|r| for the unit 


vector in the radially outward direction, this can be written 
as 


F(r) = GmM grad 


GmM 


r2 


F(r) = — e,. 

12. By definition, a spherically symmetric mass distribution 
has a mass density function p which depends on r only 
through its magnitude 


[nl] = 4 fet + oy + e, 
Of the functions listed, only (c), (d) and (f) meet this 
criterion. : 


13. From Equation (19) of Section 2, the gravitational 
potential energy of the particle of mass m is 


GmM 


This is the same potential energy as was found in Exercise 11 
for gravitational attraction by a single particle of mass M at 
the origin. Hence the gravitational force on the particle of 
mass m is, as before, 


F(r) = — grad U(r) = eels, 


p2 


e;. 
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14. Suppose that Body 1 has mass M and that Body 2 has 
mass m. Using the given hint, the motion of the centre of 
mass C’2 of Body 2 (which is also its geometric centre, by 
spherical symmetry) is the same as that of a particle of mass 
m acted upon by the same total force. However, we know 
from the result of Exercise 13 that the gravitational force 
experienced by such a particle due to Body 1 is directed from 
C2 to the centre C; of Body 1, with magnitude GmM/R?. 
This is therefore the total gravitational force exerted by 
Body 1 upon Body 2. 

By a similar argument, or by application of Newton’s third 
law, Body 2 exerts upon Body 1 a force of the same 


magnitude but in the opposite direction, that is, from Ci to 
C2. 


15. (i) The magnitude of the gravitational force on a 
particle of mass m at the Earth’s surface (where r = R) is 
mg (from the simple model) and mMG/R? (from Newton’s 
universal law of gravitation). Since these must be equal, we 
have g = MG/R’. 

(ii) From the given values, we obtain 

gh” - 98tx (G4Ax 17 


M = =e 
G 6.57 x10 = 


~ 6.0 x 10** kg. 
16. At height h above the surface, the gravitational force on 


a particle of mass m has magnitude mMG/(R + h)’, so that 
the gravitational acceleration at this height is 
MG 
CR hy) 
At ground level, we have 
MG 
wer 
from Exercise 15(i). The gravitational acceleration at height 
h is therefore 
GR: ai Sthtegs 
(R+h)? (1+h/R)?’ 
With g = 9.81ms *, h = 10*m and R=6.4 x 10° m, this 
gives the value 9.78ms ° for the gravitational acceleration 
at a height of 10 km. 


Solutions to the exercises in Section 3 


1. Forces (a), (b), (c) and (f) are central forces, since these 
are directed along r or along —r, depending on the sign of 
the constant K (and, in the case of (c), on the sign of x). 


2. From the result of Exercise 13 in Section 2, the 
gravitational force acting on the particle of mass m is 
GmM 
ri = 3) 


giving g(r) = —GmM/r°. 


3. (i) By the definition of the cross product of two 
vectors (Unit 14 Subsection 3.5), 


L =r x me 
is perpendicular to r whenever L + 0. 


(ii) Combining Newton’s second law and the definition of a 
central force (Equation (1) of Section 3), we have 


mr = F(r) = g(r)r, 


for some scalar function g(r). Then we obtain 


=rxmr+r xX mr 

=0+rxg(r)r=0, 
where we have twice used the fact that a x a = O for any 
vector a. Hence L is a constant vector. 


(iii) Referring once more to the definition of the cross 
product, if L = 0 and r 4 0,r # O, then r must be in the 
direction of r or of —r. 
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4. (i) The angular momentum of the particle is 
br % eee. 
but since the motion is in the (2, y)-plane, we may apply the 
equations 
r=re, and r=7e,+rep, 
(Equations (8) and (12) of Section 1). Since e, x eg = k, we 
obtain 


L= re, xm (Fe, + reo) 


= mrre, X e; + mre, xe, = mr? 6k. 
Hence L = Lk, where L, = mr?6. 
(ii) Since m is constant, L, = 0 is equivalent to 
d / 2: 
dt (r 0) = 0, 
which is Equation (3) of Section 3. 


5. (i) Isotropic central forces have the form 

F(r) = f(r)e-. Of the forces listed in Exercise 1, the 
following are isotropic central forces: 

(a) since F(r) = Kre;; 

(b) since F(r) = (K/r)e;; 

(f) since F(r) = Ke,. 


(ii) For the gravitational force we have, from Exercise 2 (or 
from Exercise 13 in Section 2), 


GmM GmM 
F(r) = — ee 
so this is an isotropic central force with f(r) = —-GmM/r?. 


6. The effective potential energy function is 


2 
ae Se GmM 
2 ee 2mr2 Eos 
Differentiation of this equation gives 
ge ie , GmM 
dr = — mr ee 
Now dU" /dr = 0 at r=r’", which leads to the solution 
L? 
~ m2MG 
The corresponding effective potential energy is 
: L? GmM °M?G? 
Usin = US (e*) = SE = 
2m(r*)? — pf Be 


7. If E=Uc4. then we have * = 0 and r =r" at all times. 
This orbit is therefore a circle of radius r*. 


8. (i) The angular momentum is L, = mr?6. At the 
instant of launch we have r = R and, since there is then no 
radial component of velocity, the = speed is v = R86. 
Hence L. = mvR. 


(ii) The total mechanical energy F is the sum of the kinetic 
and potential energies, and is constant. Initially, the kinetic 
energy is mv" and the potential energy is —-GmM/R.. 
Hence, for all subsequent times, 


M 
i= ime" = skied 
(iii) The motion is bound provided that E < 0. So for 
bound motion, we must have 


5 MU Se or 


2MG 
2 a R 


9. (i) Since v? = MG/(2R) the total mechanical energy 
is, from Exercise 8(ii), 
Pie Pre oe GmM = _ 38GmM 
- R AR 


Since F is negative, the motion is bound. 
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(ii) In general, bound motion has two turning points 
(although for the special case of motion in a circle the radius, 
r, is constant). From the given hint, the equation for turning 


points is E = U°#(r), where 
Sis GmM 


eff 
nee Fa ‘i 
From Exercise 8(i), we have L. = mvR, so the equation for 
Tmin and rmax is 
(mvR)  GmM_ _  3GmM 
ee ee 
Putting s = R/r, the equation becomes 
mv*s* GmMs © Sih 
ee eae ee 
The given condition v? = MG/(2R) simplifies this to 
is —s= —3, or s*—4s+3=0. 
The solutions are s = 1 and s = 3, which correspond to 
Tinax = J one Te = +R. 
(We must assume here that the spherically symmetric body 
at the origin has a radius less than +R. Otherwise the 
orbiting particle would collide with it!) 


10. (i) The equation of orbits is 
: = 1+ ecos(0@ — 6). (1) 
- 


The right-hand side of this equation takes values between 
1 —e and 1 +e (inclusive) as 6 varies. The minimum value of 
r corresponds to the maximum value of the right-hand side, 
which is 1+ e. Hence we have 
l 
lt+e 
(ii) Ife <1 then 1—e > 0, so that the right-hand side of 
Equation (1) is always positive. Its minimum value 1 — e 
corresponds to the maximum value of r, giving 
l 
1 ert: 
(iii) If e > 1 then 1 — e < 0 and the treatment of part (ii) 
cannot apply, since r is always non-negative. The condition 
r > 0 is equivalent to 
1+ €cos(? = @p)-> 6, 
which is satisfied by the range of angles 
09 — 9a <O0<609+ 62, 
where 6, = arccos(—e~*). For this range of values of @ there 
is no upper bound for the values of r, so the motion is 


unbound in this case. The corresponding orbits (for e > 1) 
are hyperbolas, as illustrated in the diagram below. 


Timm = 


Tmax = 


Figure 1 


(iv) Orbits are bound for e < 1 and unbound for e > 1. 
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11. When e = 0 we have (from Equations (18) and (16) of 
Section 3) 
L: 
m2MG 
This is a circle of radius |, for which rmax = Tmjy = !. From 
Example 1 of Section 3, the total mechanical energy is 
po ee mS | 
21 55 eee 
(The radius of the circle agrees with the value of r* found in 
Exercise 6, as expected from Exercise 7. Since 7 = 0, the 
corresponding energy is # = U daar and the value of F just 


obtained agrees with that for U®"(r*) found in Exercise 6.) 


r—-—l= 


12. (i) From Equation (13) of Section 1, the acceleration is 
i = (# — rO”)e, + (276 + rO)eg. 

Here r = d and 6 are constants, so that 7 = * = 6 =0. Hence 
#¥ = —d6’e,, 

showing that the acceleration is in the inward radial 

direction, with magnitude |r| = dé”. 

(ii) From Equations (8) and (12) of Section 1, we have 
r=rTe, and r=tTe,; + rbep = roe, 

so that the angular momentum is 
L=r x mr = mr76e, x eg = md?6k. 

Thus the magnitude of the angular momentum is L. = md?6. 


(iii) The gravitational force on the satellite is (from 
Exercise 13 of Section 2) 


F- GmM  GmM 
ee oe d? oe 
so that, from part (i), mr = F gives 
MG 
de? = a (>) 
Using the result of part (ii) to eliminate 0, this becomes 
| Fe £6: eee 
m2d3 . d? ~ m2MG’ 


in agreement with the value of the orbital radius / found in 
Exercise 11. 


(iv) The gravitational potential energy is -GmM/d. The 
kinetic energy is (from Equation (2)) 


im(d0)? = ane 
so that the sum of the kinetic and potential energies is 
GmM 
oS oe 
With d = 1, this is the same as the expression for EF found in 
Exercise 11. 


13. (i) Ifthe Earth is assumed to be spherically 
symmetric, and d is the distance between the Moon’s centre 
of mass and the Earth’s centre, then by Newton’s universal 
law of gravitation, the magnitude of the gravitational force 
exerted on the Moon is |F| = GmM/d?. 


(ii) We have |#| = dé. The angular velocity 6 is constant, 
and therefore equals 27/7, where T is the orbital period. It 
follows that 


And 

fee 

(iii) The result g = MG/R? is shown as in the solution to 
Exercise 15(i) of Section 2. 


= 


52 


(iv) From Newton’s second law mr = F applied to the 
gravitational motion of the Moon, and from parts (i) 
and (ii), we have 
4n*dm _ GmM 
To Sal 2 
which is equivalent to 


T? = (a) d°. 


From part (iii), this is the same as 


oa. 4n? 3 
= (ia)¢ 


as required. 


(v) Using d= 60R, g =9.8ms *, T = 28 x 24 x 3600s and 
R=5.3 x 10° m gives 

dad An X OO XE BK 

gR?T2 9.8 x 282 x 24?.x 36002 5 — 
This is the ratio of the two sides of the equation in part (iv), 
which should be close to 1 if the prediction is accurate. (If 
R= 6.4 x 10° m is used instead of 5.3 x 10° m, with the other 
values unchanged, then the value of this ratio becomes 0.95.) 


14. (i) Ina time 7 the whole area of the ellipse is swept 
out. Since A is constant, we then have 


Squaring both sides of this equation, and putting 


b> =al= aL: 


we obtain 


as required. 


(ii) For the circular orbit of Exercise 13, the acbital radius 
isa =b=d. Use of MG = gR? (from Exercise 13(iii)) then 


gives the result of Exercise 13(iv) once more. 


Solutions to the exercises in Section 5 


1. (i) Ifthe orbit is bound, then we have 

mMG 
21 ( 

together with 


C= 


= 4) 


b=:— gnd <= 
a 


(Equations (4) and (5) of Section 2). Hence 


e? —1 = —b?/a? = —I/a, giving 
EE _mMG 
2a 


(ii) For any orbit, Exercise 10(i) of Section 3 gives the 
distance of closest approach as 


l 


ie Se 
Hence we obtain 


f= me le — 1). 


min 
(As a check, both of these results agree with that of 
Exercise 12(iv) of Section 3 when the orbit is circular, with 
radius d = a = rpin and eccentricity e = 0.) 


min = 


MS1T204 30 Solutions 


2. The result 
er = gh’ T* 

Ar? 
can be obtained just as in parts (i)—(iv) of the solution to 
Exercise 13 in Section 3, though here d is the radius of the 
satellite’s orbit rather than that of the Moon’s. On putting 
g=9.81ms 7, R= 6.4 x 10° m and T = 1day = 24 x 3600s, 
this gives 


ao 


1 
Ar heed 
or d~ 4.24 x 10'm~6.6R. 


The corresponding value for the height above the Earth’s 
surface is 


d — R~ 3.6 x.10' m. 
(The use of geostationary orbits for communications 
satellites was first suggested by the science fiction author 


Arthur C. Clarke.) 


(9.81 x (6.4 x 10°)” x (24 x 3600)”) 


3. (i) For bound motion, we have both 


b= ime" — “= 
and (from Exercise 1(i) of this section) 
a _mMG 
2a 


On equating these expressions and solving for a, we obtain 


= 2 vy aa 
~ th. Mo 


(ii) From oo ) and ( 


But also, it is oe from pera (16) of Section 3 that 
L2 
~ m2MG" 
Since L, = mvR is given in the question, we have 
v* R? 
MG 
and, using the result of part (i), 


wh? (2 v" ks 
eS oe Cee 
: | MG (F wa) 
R 


1 
v°R v? ae 
ae Se. (ie) : (sa) 


v?R 
MG’ 
where the magnitude has been taken to ensure that e > 0. 


(Note that the condition e < 1 must correspond to the given 
condition v? <2MG/R for bound motion.) 


Alternatively, this result could have been obtained from the 
formula 


) of Section 2, we have 


= 


; os 


or € = 


GmM 


ie = 51 fe. 
(iii) For circular motion, we have e = 0 and 
2. MG 
p 
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